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Abstract 

We prove regularity for a class of boundary value problems for first 
order elliptic systems, with boundary conditions determined by spectral 
decompositions, under coefficient differentiability conditions weaker than 
previously known. We establish Fredholm properties for Dirac-type equa- 
tions with these boundary conditions. Our results include sharp solvability 
criteria, over both compact and non-compact manifolds; weighted Poincare 
and Schrodinger-Lichnerowicz inequalities provide asymptotic control in 
the non-compact case. One application yields existence of solutions for 
the Witten equation with a spectral boundary condition used by Herzlich 
in his proof of a geometric lower bound for the ADM mass of asymptoti- 
cally flat 3-manifolds. 
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1 Introduction 

Elliptic systems based on the Dirac equation arise frequently in problems in 
geometry and analysis. Applications to positive mass and related conjectures 
in general relativity motivate this paper, and involve boundary value problems 
on compact and non-compact domains [31,37,38,64]. 

Previous existence and regularity results [6, 12, 16, 40, 50] are insufficient 
for these applications, for various reasons. The Agmon-Douglas-Nircnbcrg ap- 
proach based on freezing coefficients and explicit kernels for the constant coeffi- 
cient inverse operator, leads only to boundary conditions of Lopatinski-Shapiro 
type [44]. The pseudo-differential operator approach [40,56] handles non-local 
boundary conditions such as the spectral projection condition of Atiyah-Patodi- 
Singer [6] , but the assumptions of smooth coefficients and product-type bound- 
ary metric [6, 12, 16, 40] are unnatural and, as we shall show, unnecessary. 

In this paper we provide an essentially elementary proof of existence and 
regularity for first order elliptic systems with "Dirac-type" boundary value con- 
ditions. These encompass both pointwise (Lopatinski-Shapiro) and non-local 
(spectral) boundary conditions, and do not require product metric structures 
on the boundary. We obtain explicit necessary and sufficient conditions which 
ensure the solvability of natural inhomogeneous boundary value problems, over 
both compact and non-compact manifolds with compact boundary. 

The coefficient regularity conditions, for both the elliptic system and the 
boundary conditions, are rather general. For example, they are weaker than 
those in the pseudo-differential operator approach of Marschall [48] . It seems 
likely that the boundary conditions will admit some generalizations; the bound- 
ary data is H^/'^ whereas there are recent results for a certain constant coefficient 
Dirac equation with boundary values on a Lipschitz hypersurface [1] . 

Note that there is an extensive literature on applications of Dirac operators 
to index problems on compact and non-compact manifolds [12, 18] which we do 
not address, although many aspects of our results are no doubt relevant to such 
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applications; the results here are focussed on applications to energy theorems 
in general relativity. 

The motivating example of the Dirac (Atiyah-Singer) operator is described 
in some detail in §2, where the Schrodinger-Lichnerowicz identity with suitable 
boundary conditions combines with a Lax-Milgram argument to reduce the 
existence question to that of showing that a weak (L^) solution of an adjoint 
problem is in fact a strong {H^) solution. This weak-strong regularity property 
turns out to be the key technical step, and the focus of much of the paper. The 
difficult case is regularity at the boundary; interior regularity is established in 
§3 using standard Fourier techniques, for general first order elliptic systems. 

§4 reviews conditions under which a symmetric operator has a complete 
set of eigenfunctions; these are used to to control the boundary operator in 
later sections. In §5 we prove regularity results at the boundary, for a class of 
operators much broader than Dirac equations, with weak assumptions on the 
continuity /regularity of the operator coefficients. The main technical tools are 
the identity (5.15), and some basic spectral theory. The boundary conditions 
of §5 follow from the requirements of the arguments of the regularity theorem, 
and some additional work is required to apply them to first order systems. This 
is carried out in §6, for equations of Dirac-type near the boundary, for which 
the boundary operator is self-adjoint. The resulting boundary conditions are 
naturally presented in terms of graphs over the space of negative eigenfunctions 
of the boundary operator. 

The boundary value problems considered have a Fredholm property, and 
admit an explicit solvability criteria involving solutions of the homogeneous 
adjoint problem. These properties are established for compact manifolds with 
boundary in §7, and for a large class of non-compact manifolds with boundary 
in §8. The analysis of the non-compact case relies on two a priori inequalities: a 
weighted Poincare inequality, and a Schrodinger-Lichnerowicz inequality. These 
inequalities imply the manifold is non-parabolic at infinity in the sense of [18]. 
The weighted Poincare inequality is established in §9 in a number of cases, in- 
cluding the important cases of manifolds with asymptotically flat or hyperbolic 
ends. The Schrodinger-Lichnerowicz inequality follows in applications from an 

estimate derived from an identity of Schrodinger-Lichnerowicz type. 

In section 10 we show that common pointwise and spectral boundary con- 
ditions for the Dirac equation are elliptic in the sense of our conditions. These 
calculations form the basis for §11, which verifies several positive mass theo- 
rems [31,37,64]. Appendix A collects some relevant properties of tensor and 
spinor fields on manifolds with W'^'^^'^ differentiable structure and W'^''^ metric, 
k > n/p. 

2 The model problem 

In this section we use the Riemannian Dirac equation to illustrate and motivate 
the existence and regularity results of the following sections. 

Consider an oriented manifold M with Riemannian metric g and a represen- 
tation c : Ce{TM) End(S') of the Clifford algebra Ct{TM) on some bundle 
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S; with our conventions, 

c{v)c{w) + c{w)c{v) = —2g{v,w) . 

Clifford representations are discussed in detail in [5,45]. S carries an invariant 
inner product, {c{v)ip , c{v)-ip) = = I^P IV'P) with respect to which 

c{v) is skew-symmetric, for all vectors v. 

A Dirac connection [12, 45] is a connection on the space of sections of S 
which satisfies the compatibility relation 

c?(</>, c{v)^) = c{v)^l;) + {cl>, c{v)Vij) + (</., c{Vv)ij) , (2.1) 

where V also denotes the Levi-Civita connection on vector fields. 

Spin manifolds provide the fundamental example, with S a bundle of spinors 
associated with a Spin principal bundle which double covers the Riemannian 
orthonormal frame bundle. In this case there is a covariant derivative V de- 
fined in terms of a local orthonormal frame ejfc, k = 1, . . . ,n, with Riemannian 
connection matrix uiij{ek) = g{ei,Vei,ej), by 

Ve, V = De,^' <t>I - W , (2.2) 

where if; = ijj^ and <^i, I = 1, . . . ,dimS', is a choice of spin frame associated 
with the orthonormal frame e^. The expression (2.2) may be abbreviated to 
V = d — ^uJijC^eK Note that there are other examples of Dirac bundles and 
connections, eg. [45, example II.5.8]. 

The Dirac operator of a Dirac connection V is 

PV^ = c(e*)Ve,^ ; (2.3) 

in the spin case this is sometimes called the Atiy ah- Singer operator. When the 
spinor representation is irreducible'^ , a classical and very important computation 
[55] shows that 

= v*V^ + \R{g)^ , (2.4) 

where R{g) is the (Ricci) scalar curvature of g. This leads to the Schrodinger- 
Lichnerowicz identity [46, 55] 

(|V^|2 + lR{gM\' - \V^pf) *l = d ((V, {ciaej) + 5^,-)V^'V) *e') , (2.5) 

which when integrated over the compact manifold M with boundary^ Y be- 
comes 

/ (I VV;|^ + lR{gMf - = / c(ne^)VAV) • (2-6) 

JM JY 

Here n is the outer normal vector at y = dM and {e^} is a compatible or- 
thonormal frame on Y. The boundary term may be simplified by introducing 
the boundary covariant derivative 

V = d - \LOABc{e^e^) , 

^Reducible representations lead to interesting formulas with ^R{g) replaced by more com- 
plicated curvature endomorphisms, c.f. §11 

^Throughout this paper we use the geometer's convention, that a manifold with boundary 
contains its boundary as a point set. 
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and the boundary Dirac operator^ 



VyiJ = c(ne^)VAV' • (2.7) 
Denoting the mean curvature hy H = Hy = g{n, Vg^e"^) gives 

c{ne^)VA4>) = 'Dy4> + \Hij) . (2.8) 

We use conventions which give H = 2/r > for M = M.^ — B(0,r), the exterior 
of a ball of radius r, with the outer normal n = —dr- If a; is a Gaussian boundary 
coordinate (a; > in M, x = on y and dx = —n), then near the boundary we 
have 

Vi/j = -c{n){dx + Vy + \H)iI) . (2.9) 

We now seek boundary conditions for which the equation Vif) = / is solvable, 
following a well-known argument [32, 37, 51]. Suppose M is a compact^ manifold 
with non- negative scalar curvature, R{g) > 0, and /C : H^^'^iY) H^^'^iY) is 
a bounded linear operator such that 

j> (-0, Vy^) + iHil)) < whenever Ki/; = . (2.10) 

Suppose further that M admits no parallel spinors. Define the space H^{M) 
as the completion of the smooth spinor fields with compact support (in MUY) 
which satisfy the boundary condition ICtp = 0, in the norm 

MlliM) ■■= Ij'^^l' + iRiaMl") ■ (2-11) 

The boundary condition (2.10) combined with the Lichnerowicz identity (2.6) 
and the curvature condition R{g) > now ensures that the bilinear form 

Jm 



IM 

is strictly coercive, a{il),'4)) > Wi^W^H^ (M)' ^'^^ spinor field / G L'^{M), the 
linear functional i— s- Jj^{f,'D(j)) is bounded on Hi^{M). Coercivity and the 
Lax-Milgram lemma show there is a unique ijj G H^{M) such that 

Jm 

and we would like to deduce that = /■ Now * := Vijj — / G L'^{M) is a 
weak solution of the Dirac equation; that is, 



/ 

Jm 



V(f)) = V (/) G H]c{M) . (2.12) 
i_M 

^Both — > c(e^) and c{ne^) give representations of the Clifford algebra of the 

boundary tangent space; the choice of c{ne^) is made hero for convenience [37]. 

*The case of an [M, g) which is asymptotically flat with compact interior, which is relevant 
to the positive mass theorem, is discussed along similar lines in §8 and 11. 
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If we could show that * is in fact a strong solution, that is, * G H^{M), then 
we could integrate by parts to conclude 

/ (P*, </.)+/ c{n)cj)) =0 Vc/. G HjciM) , 
Jy 

and thus = and ^(*, c(n)0) = for all 4> G iJ^/^^y) g^^^.^^ ^j^^^^ ^^^^ ^ q. 
This would give the boundary condition ^'|y G c(n) (ker /C)-*", which we suppose 
may be re-expressed as /C^* = 0, for some "adjoint" boundary operator /C. This 
would give ^' G so if finally we suppose that /C also satisfies the boundary 
positivity condition (2.10), then we could conclude from a(*,^') = and the 
coercivity of a(-, ■) with respect to the norm || ■ \\ffi , that * = as desired. 

K 

The key technical difficulty in this classical argument lies in establishing the 
"Weak-Strong" property, that weak (L^) solutions lie in H^. In the following 
sections we will prove this property for a large class of elliptic systems, under 
rather general boundary conditions; see §5 and §6. 

Two model boundary operators illustrate the possibilities for achieving the 
required conditions. The APS (or spectral projection [6]) condition arose in 
Herzlich's work [37]: 

/C = P+ , (2.13) 

where P+ is the (M)-orthogonal projection onto the positive spectrum eigenspace 
of the boundary Dirac operator Dy Using the relation c{n)'DY = —'DYc{n), 
which shows that the spectrum of Dy is symmetric about G M, we find that 
/C = /C, provided there are no zero eigenvalues. 

The eigenvalue estimate for F 2± 5^ of Hijazi and Bar [7, 39] 



\X(Py)\ > V4vr/area(y), (2.14) 

shows that in this case there are no zero eigenvalues. In addition, if we have 
the mean curvature condition 

Hy < Vl67r/area(y), (2.15) 

then /C (and /C) will satisfy the boundary positivity condition (2.10). In con- 
clusion, if y = dM ~ S"^ satisfies (2.15), then (assuming the Weak-Strong 
property can be established) the above argument shows ^' = and thus the 
equation Vijj = f with boundary condition P+ip = is uniquely solvable, for 
any / G L''{M). 

The chirality condition was used in [30, 31]. For a slightly simplified version 
of [31], suppose M is a totally geodesic hypersurface in a Lorcntz spacetime, 
with future unit normal vector cq, and consider the connection on spacetime 
spinors, restricted to M. Along Y = dM we define 

e = c(eon), (2.16) 

which satisfies the chiral conditions 



= 1 , ec(n) -|- c(n)e = , 

(<^,e^) = (e(/),V), eVy + Vye = 0, 



(2.17) 
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and then the boundary operators 

/C± = i(l±e). (2.18) 

Assuming either of the two conditions IC±ip = gives etp = which imphes 

{tp,VYtp) = =f{ij,VYeiP) = ±{ij,eVYip) 

= ±{e^,VY^) = -{^,Vy^)=0. 

If wc further assume that Hy < then (2.10) follows directly. In general 
relativity the condition Hy < is the defining property for y to be a trapped 
surface. 

Since the /C-t's are complementary orthogonal projections, we have (kerlC^)^ = 
ker/Ciji, so V £ c(n) (ker/Ci)-*- exactly when c{n)ijj G ker/Cip, which gives 
■0 G ker/C-j-, and /C± = IC±. In this case we conclude (still assuming the Weak- 
Strong property can be established) that if Hy < then Vip = f, with either 
of the boundary conditions eijj = ±-0, is uniquely solvable. These examples are 
discussed further in SIO. 



3 Interior Regularity 

In this section we establish regularity away from the boundary for weak (L^) 
solutions of first order elliptic systems. We consider equations of the form 

Cu := a^djU + bu = f , (3.1) 

where u, f are sections respectively of A'^-dimensional real vector bundles E, F, 
both over an n-dimensional manifold M without boundary^, and a^,b, j = 
1, . . . , n, are sections of the bundle of endomorphisms of £^ to F. We assume 
that E,F are equipped with fixed smooth inner products, denoted by (•,•). 
The length determined by (•, •) will be denoted invariably by \u\'^ = {u,u). To 
simplify notation, the respective bundles usually will be understood, and thus 
L'^{M) will generally mean L'^V{E), the space of sections of E, or L^r(F), 
depending on context. 

Remark 3.1 There is no loss of generality in considering real bundles, since 
complex and quaternionic bundles may be viewed simply as real bundles with 
additional algebraic structure. For example, a Hermitean vector space of di- 
mension n is equivalent to a real vector space of dimension 2n with a skew 
endomorphism J satisfying = —1, with the Hermitean inner product ( , ) 
and real inner product ( , ) related by (u, v) = {u, v) — i{u, Jv). 

Define the indices 2 = 2(n), n* = n*(n) by 

2 = n* = n for n > 3, 

2 = 106, n* = ^ for n = 2, (3.2) 

2 = GO , n* = 2 for n = 1 , 

^If M has boundary dM ^ 0, then the interior M = M — dM is a (noncompact) manifold 
without boundary, to which the results of this section will apply. 
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where 10^ represents any large constant. Note that if M admits a Sobolev 
inequahty with constant Cs 



u\\j^2 < C5||u||^i, 



then we also have 



\\fu\\L2<Cs\\f\\L.'\\u\\m ■ 



(3.3) 



Another basic fact is the inequality 



WfgWw^^r < c 



which shows that VF^'" n (in particular) forms a ring under addition and 
multiplication of functions. For 77, = 1,2 the C° is superfluous here, of course. 

With one exception, it suffices to assume throughout that the underlying 
manifold has a C°° diffcrcntiablc structure. The exceptional point arises in §10 
in the construction of approximately Gaussian coordinates in a neighbourhood 
of the boundary, when the metric has low regularity. The description in Ap- 
pendix A of W'^'^^'P differential structures, k > n/p, establishes the necessary 
consistency conditions in this case. 

We will assume that , b satisfy the regularity conditions 



The conditions (3.4) are preserved by bundle frame changes in W^'"'* fl C°, 
by the above ring property. In particular, even if the bundle metrics (•, •) on 
E,F are only in W^'^''* Pi C^. by the Gram-Schmidt process wc may construct 
W^''"* n frame changes which make the metric coefficients constant. Since 
this changes the operator coefficients , b respectively by W^'"'* DC'^ , W^* affine 
linear transformations, there is no loss of generality in assuming the metrics on 
E, F to be locally constant. 

The conditions (3.4) mean that M can be covered by open neighbourhoods 
ffa with VF^'"* n C° bundle transition functions, such that the local coefficients 
a^, b satisfy the stated regularity. Frame changes satisfying Sobolev conditions 
are also discussed in detail in Appendix A. 

We require that satisfy the ellipticity condition, that for each p & M 
there is a coordinate neighbourhood p e U C M and a constant 77 > such 
that 



for all X G U, ^ G T*M and V ^ E^, where is measured by a fixed 
background metric g, which we may assume to be C°°. Note that (3.5) implies 
the fibres of E, F must be of the same dimension. 
A weak solution of (3.1) is u G Lf^^{E) such that 



e H^/jf (M)nCO(M) , 
b e L-;(M). 



(3.4) 



(3.5) 




(3.6) 
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for all (j) € C^{M), where dvM = "jdx, 7 > 0, is a coordinate-invariant volume 
measure on M with 7 £ W^'^ (U) H C^{U) and dx is coordinate Lebesgue 
measure, in any local coordinate neighbourhood U. Here the formal adjoint 
is defined with respect to dvM and the inner products on E, F. Thus in local 
coordinates, 

/:t^ = -ta^dj(j) + Cb - j-^djCa^-f))cl) , (3.7) 

where the transposes ^a^ are defined with respect to the local framing forms of 
the inner products of E, F. 

The proof proceeds by establishing various special cases, starting with a 
constant coefHcient operator acting on sections of a trivial bundle E over the 
torus T". This type of argument is very standard. 

Proposition 3.2 Suppose u G L^(T"') is a weak solution of Cqu = f where 
f G L^(T"). £0 = o,Qdj with Oq constant and satisfying the ellipticity condition 
(3.5). Then u e ifi(T"). 

Proof: Wc regard T" = IR"/Z". Fix a mollificr <j)^ = e-''(p{{x-y)/e) G C~(M'') 
with (p{-x) = (j){x), and set u, = (j),*ue ^^(T"). Then 

f 9 

^QUe = al-^(j)e{x - y)u{y) dy 
and thus the definition of weak solution gives 

/ {tp,CQUe)dx = 4>e{x - y){CQtp{x),u{y)) dydx 

I I {-ao^My-x)ip{x),u{y))dydx 

Jjn Jjn. OyJ 

{^Qi(pe*ip){y),u{y))dy 

/T" 

{4>e*'^P{y),f{y))dy 
{i^{y),(t>e * f{y))dy . 



j 

j 
I 



2 



Thus Cqu^ = fe = (j)e * f, and we note that fe ^ f strongly in L . Now 
the ellipticity condition (3.5) and the Plancherel theorem ensure that for all 

V e /fi(T"), 



Jfn Jjn 

= r/"^ / |>Co^^|^ dx , 



and thus 



/ \due\^dx < r/-^ / \fe\^dx . 

Jin Jin 

Since u, f e , it follows that ^ u strongly in H^. 
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Proposition 3.3 Under the conditions of Proposition 3.2, the map Cq + \ : 
iJ^(T") — >■ L^(T") where A = ttt/, is uniquely invertible, and for allu G H^{T^), 



\u\ 



1^,1 < V5/r?||(£o + A)u|U2 . (3.8) 



Proof: Write u = Y^k^jn iifcC^"'^'^, where the coefficients Uk = Jj„ u{x)e~'^'^^^'^ dx 
are valued in C''^, the complexification of the real vector space modelling the 
fibres of E. We then have 

/ \{Co + A)tip = V \{27rikjai + A)ttfcP, 

and using the vector length inequality \a + fep > xl*^!^ ~ x/{^ ~ x)I^P with 
X = ^, we find 

\{2TrikjaQ + A)^^^^ > ^\27rikjaQUk\'^ — X'^\uk\'^. 

For k ^ this is greater than 7r^r7^|«fcp, whilst for A; = we have \{2mkjaQ + 
= '^^kitP = Tr^^^l'WfePj hence 

/ \{Co + A)up > V TrVl^tfcP = ttV / I^^P dx , 

which shows Co + X has trivial kernel. Choosing x = 1 — l/(2|fe|) shows in fact 
that 

\{2Trikjai + A)njtp > Tr^ri'^{2\k\ - lf\uk\'^, 

for all A; G and all Uk G C^. Since 2\k\ - 1 > 1 for all A; G Z", we obtain 
(3.8). Moreover, this shows also that the N x N complex matrices 2TrikjaQ + A 
are invertible for any A; G Z", which gives a direct construction of the inverse 
of the operator Cq + X. ■ 

Theorem 3.4 Suppose u G L^(T") is a weak solution of 

Cqu + Bou + Biu = f (3.9) 

where f £ L'^ and Cq = a^dj is a constant coefficient first order operator 
satisfying the conditions of Proposition 3.3 with ellipticity constant rj, where 
Bi : ^ is bounded, and where Bq : is a linear map satisfying 

\\B4h^^l^ < »7/3, WbIWh^^l^ < r//3 , (3.10) 

where bI is the L^iT"")- adjoint of Bq. Then u G H^{T') is a strong solution 
of (3.9), and there is a constant C, depending only on rj and ||-Bi||^2_>i2, such 
that 

II^^IIhI <C(||/|U2 + 11^.11^2) . (3.11) 
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Proof: Construct the iteration sequence w^'^'' G H^, A; = 0, 1, ... by defining 
to be the solution of 

{Co + A)u;('=+i) = -Bow^''^ + f , (3.12) 

with w^^^ = 0, where / = f + \u — Biu € by the assumptions. This 
equation with A = Trry is uniquely solvable by Proposition 3.3. The difference 
^(fc+i) ^ y^{k+i)_^{k) satisfies {Cq + \)v^^+'^'^ = -Bqv'-'') and the estimate (3.8) 
shows that 

The iteration is thus a contraction and converges in ^ to w <^ satisfying 
(£o + Bq + X)w = /, and then v = u — w ^ L? \s & weak solution of (£o + + 
X)v = 0. Now £q is also elliptic with the same ellipticity constant rj, so there is 
z G satisfying (£q + Bq + X)z = v. Since f is a weak solution, 

{iCl + Bl + X)(l>,v)dx = V<^GiJ^(T"), 



/ 



we may test with (p = z to see that J \v\'^ = and v = 0. Thus u = v + w = 
w G as required. By Proposition 3.3 and (3.10), we have 

^ll^ll^i < \\{jro + Bo + Bi)u\\L2 + \\Bou\\L2 + \\Biu\\L2 + \\Xu\\L2 

Tj 

< II/IIl2 + gll^^llf/i + (II^i||l2^l2 +??7r)||u||i2 . 

Since \/5 < 3, the estimate (3.11) follows. ■ 

Next we consider operators with non-constant coefficients. Let Cs be the 
T" Sobolev constant 

Iklli2(xn) ^ C's||'"l|i/i(T") , (3.13) 
where 2 is defined in (3.2). 

Proposition 3.5 Suppose u G L^(T") is a weak solution of the system of equa- 
tions 

Lu := a?dju + bu = f (3.14) 
where f & L?' and the coefficients G VF^'"* n b G L"* satisfy 

11^^"' 11^"* ^ 1^ ^'-''^ 

where a^, j = 1, . . . ,n, are constant matrices with ellipticity constant rj. Then 
u G (T") is a strong solution of (3.14) ■ 

Proof: It will suffice to show that C admits a decomposition satisfying the 
conditions of Theorem 3.4. Since L°° is dense in L" , for any e > we may 
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find 60 € L"*, 61 G such that 6 = 60 + ^'i and ||6o||Ln* < e- We choose 
e = r]/{10Cs)- Then Bqu := bou + (a-' — a-^jdjU satisfies 

||-Bo^x||l2 < ||6o||in* \\u\\j^2 + \\a^ - al\\L'^\\du\\L2. 
Using (3.16) and the Sobolev inequaUty (3.13) gives 

\\Bou\\l, < ^{\\u\\hi + \\du\\L2) < "^MIhi, 

so ||i3o||Hi^_L2 < ry/3. Clearly Biu := 600^ is bounded on L^, and it remains to 
verify the bound on the adjoint operator 



bIw := %w - Ca^ - ^al)djw - dj^a^ - *a^) 



w. 



Again using the Sobolev inequality and the conditions (3. 15), (3. 16) we find 

ll^dtylUz < C5||6o||in* Iklki + -^al\\L°°\\9w\\L^ + Cs\\dja'\\Ln*\\M\H^ 
^ ^11 II 

SO the conditions of Theorem 3.4 are met and the result follows. ■ 

On a general compact manifold we define the Sobolev space (M) by the 
norm 

\Mhhm)= I (|Vu|2 + |u|2)di;M , (3.17) 

where the lengths \u\^, \Vu\^ are measured using the metric ( , ) on sections of 
E and a fixed smooth background metric g on TM, and where V is a (covariant) 
derivative defined in local coordinates on M and a local framing on E by 

Vi = di-Vi. (3.18) 

We assume the charts on E,M are such that 

r, e Lf„; . (3.19) 

Note we do not require that V be compatible with the metric on E. If M 
is compact then the space H^{M) is independent of the choice of covariant 
derivative: 

Lemma 3.6 Suppose M is compact and V, V are covariant derivatives satisfying 
(3.19). Then there is C > such that for all u G H^{M), 

C-^ I (|V^x|2 + \uf)dvM < [ i\^uf + \u\'^)dvM <C [ {\Vu\'^ + \uf)dvM ■ 

JM JM JM 

(3.20) 

Moreover, there is a constant Cg, depending on M, V, such that 

/ \u^dvM <Cs {\Vuf + \uf)dvM . (3.21) 

JM ) JM 
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Proof: There is a finite covering of M by charts Ua with a corresponding 
partition of unity Using the Sobolev inequahty for Ua C M", in each chart 
we may estimate the locahsation Ua = (paU by 

/ \dUafdx<C (|VUap + |naP)dvM , 

where C depends also on the decomposition T = T°° + F"* G L°° + L""* , with 
r"* smah. Again using the M" Sobolev inequality and F, F G L"* we have 

/ {\Vu\^ + \u\'^)dvM < 

JM 

< 
< 

from which the equivalence of the norms follows easily. The Sobolev inequality 
follows from very similar arguments. ■ 

We may now complete the proof of interior regularity. 

Theorem 3.7 Suppose M is a C°° n- dimensional manifold without boundary, 
and E,F are real vector bundles over M, each with fibres modelled on M'^. 
Suppose u G Lf^^{M) is a weak solution of Cu = f , where C is a first order 
operator satisfying the conditions (3.4,3.5). Then u G H^^^{M) and u is a 
strong solution of Cu = f . Moreover, if M is compact there is a constant 
C > 0, depending on , b and T, such that for all u G H^{M), 

M\m{M) < C{\\£u\\l2(^m) + MIl^m)) ■ (3.22) 

Proof: Since jC is locally of the form Cu = a^dju + bu with G W^^f H 
b G Lf^^, for each p e M there is a coordinate neighbourhood U and a constant 
rj > such that r] is the ellipticity constant of Cq = (p) , and with respect to 
the local trivialisation of — ?7 x we have the bounds 

Tj 

\y-ao\\L'^{u) < 

l|5i«'llL"*(C/) < 

where we assume without loss of generality that U = Qr = (0, i?)"" is a cube 
of side length R < 1. By paracompactness there is a locally finite countable 
covering {pa, Ua}aei. of M by such charts, with a subordinate partition of 
unity {0a}aez- Noting that supp {(pau) Qr and that (j)aU satisfies 

C{(l)aU) = (paf + a^dj(paU 

weakly, we see that it suffices to consider the case where supp u (s Qr. Assum- 
ing this, rescaling by y = x/R,x G Qr and defining ■u(y) = u{x), f{y) = Rf{x), 



cY] / (|v«„p + Kp)dx 
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a^{y) = a^{x) and b{y) = Rb{x), it follows that u G L^(T") is a weak solution 
of 

d 

^•'^"(y) + KvMy) = f{y)- 

In particular we have b G L"* (T") and 

\\d^ - a^llioo(Trn) = \\a-' - a^||ioc(Q^) < ry/lO, 

The conditions of Proposition 3.5 are satisfied, so u G ii"^(T") and thus u G 
HUM)- 

When M is compact there is a finite covering by charts {pa,Ua}, and by 
Theorem 3.4, in each chart we may estimate the localisation Ua = 4'aU by 



a\\H'^{Uo,) - Ca(||'C'"a||L2(M) + l^a ||l2(M))) 



where the H^{Ua) norm uses the coordinate partial derivatives di in Ua- To 
estimate Jjj \Vua\^, note that the Sobolev inequality (3.3) in Ua gives 

/ \TUa\'^dvM < C'l|r||^„*([;^) / {\dUa\'^ + \Ua\'^) dvM , 

so by the H^{Ua) estimate we have 

||'"a||//i(M) < CaiWCUaWh'^^M) + \\'U'a\\L'^{M)), 

for some constant Ca depending also on ||r||2,n*((7^). Since n = and 
£mq = (paCu + dj{(l>a)a-' u, with |9(/>a| < c and |(/>a| < 1, the estimate (3.22) 
follows easily. ■ 

The constant C of (3.22) can be controlled by \\a^ Ww^'P^ W^Wlp for any > n*, 
or by otherwise controlling the decompositions dja^,b G L°° + L"*. 

Higher regularity follows easily from Theorem 3.7 by a standard bootstrap 
argument: 

Theorem 3.8 Suppose u G L'^^^ is a weak solution of Cu = f in the situation of 
Theorem 3. 7, where the coefficients of Cu = f satisfy the regularity conditions 

a3 e W.'f n CO, be * , and f e , (3.23) 

for some integer k > 1. Then u G H^^^. If M is a compact manifold without 
boundary then there is a constant C = C{k,C), depending on k and \\a^\\\Yk,n* , 
Wb^Wwi^'^*' such that 

lklli?'=+i(M) < C!{\\f\\Hk(^M) + IHlni'iM)) ■ (3-24) 

Thus for any u G L^{M) such that Cu (defined weakly) satisfies Cu G H^{M), 
we have 

\H\Hf'+^{M) < C{\\Cu\\Hk{^M) + II^^I|l2(m)) • (3.25) 
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Proof: For simplicity we first treat the case k = 1. Theorem 3.7 shows u G 
^loc ^° vector of first derivatives du G L'^^^ itself is a weak solution of the 
system of equations 

Cdu + d{a^)djU = df -d{b)u . (3.26) 

Since da^ G L"* and 

\\dib)u\\L2 < C\\b\\wi,n*\\u\\Hi, 

so ||5(6) it||j^2 is bounded, this system satisfies the conditions of Theorem 3.7, 
hence u G Hf^^. The general induction step applies a similar argument: if the 
result is established M k < K — 1, and if Cu = f with coefficient conditions 
(3.23) with k = K, then du satisfies an elliptic system (3.26) of the same form 
with coefficient conditions (3.23) with A; = i^ — 1, so by induction du G H^{M) 
and thus u G H^^^{M) as required. The estimates (3.24), (3.25) follow easily 
by a similar argument and Theorem 3.7. ■ 

The coefficient conditions in Theorem 3.8 are not optimal in most cases. 
For example, if n = 3 then b G W^'^ suffices to show u G (rather than 
b G W^'^). This follows by interpolation, 

||56u||i2 < ||56||i2 ||u||loo 

< e\\u\\H2 + C(e, ||56||i2)||u||i2 , 

which shows that dbu may be thought of as the sum of a small second order 

operator, and a large bounded operator on L^. The small operator term may 
be absorbed as a perturbation of C, and the remainder contributes to the right 
hand side source term. 

4 Spectral Condition 

In this section we review conditions under which an operator will have a com- 
plete set of eigenfunctions. These conditions will be used in §5 to analyse 
boundary conditions, and thus the case of most interest concerns operators on 
a compact manifold without boundary, and in particular the first order elliptic 
systems considered in §3. However, the main result. Theorem 4.1, is stated 
in slightly more generality, which could be used to extend the eigenfunction 
representation to operators on manifolds with boundary. 

Let H he a closed subspace of W^''^(Y), with the induced norm, where Y is 
a compact manifold perhaps with boundary, and as in §3, it is understood that 
these spaces refer to sections of a (real) vector bundle E over Y. 

The abstract spectral theorem for the map A: H ^ ^'^(X) uses the follow- 
ing conditions: 

(CO) A: H ^ ^'^(X) is linear and bounded in the W^''^ topology on H. 

(CI) The Carding inequality holds: there exists a constant C such that for all 
ip & H we have 

m\% < C j^{{A^I;,A^) + {^l;,^))dvY . (4.1) 
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(C2) Weak solutions are strong solutions ("elliptic regularity") 
satisfies 

J {Aip,(j))dvY = , yipeH, 

then (pe H. 
(C3) A is symmetric: 

y (f),tpeH J {A(f), ^) dvY = J {(t>, Ai;) dvY . (4.3) 

(C4) density: 

H is dense in L'^{Y). (4.4) 

Note that in the case dY ^ 0, the space H must incorporate boundary condi- 
tions, and these will play an important role in verifying (C2), as will be seen in 
§5. 

The main result of this section is the following: 

Theorem 4.1 Under the conditions (C0)-(C4), there exists a countable or- 
thonormal basis of consisting of eigenfunctions of A, with eigenvalues all 
real and having no accumulation point in R. 

Proof: Let Ker{A) C H he the kernel of A; it is a standard fact that Ker(^) 
is finite dimensional when the Garding inequality holds — we give the proof 
for completeness. Let / < oo, be an L^-orthonormal basis of Ker(^), 

the equation Aipi = together with (4.1) shows that {■i/'i}|=i is bounded in 
W^''^. The Rellich theorem [35, Theorem 7.22] implies that from the sequence 
-02 wc can extract a subsequence ijjij converging strongly in L^, weakly in W^''^. 
The Garding inequality (4.1) with ^Jj replaced by — shows that ipi. is 
Cauchy in Vl^^'^, hence converges in norm to some G W^''^. By continuity of 
A, condition (CO), we have Aif; = 0, by continuity of norm on W^''^ it holds 
that ||i/'||l2 = 1, and it easily follows that ip € {ipill^i- have thus shown 
that {V'i}f=i is compact, which yields / < oo, as desired. 
Let now 

H = {tlj EH ■.W(t)e Ker(A) J (0, V') dvy = 0} . 

For (j) £ L'^ the map H 3 'tp ^ Jy{(p,ip) dvy G M is continuous in the 
topology (and therefore also in the W^''^ topology), thus H is closed (being an 
intersection of closed spaces) , and hence a Banach space. We note the following: 

Lemma 4.2 There exists a constant C such that 

y^PeH < CP0IU2 . (4.5) 

Proof: Suppose that this is not the case, then there exists a sequence ijjn & H 
such that 

ll^'nlU^ > nWAiPrih^ ■ (4.6) 



: If G L^{Y) 
(4.2) 



16 



Rescaling tpn if necessary we can without loss of generality assume that 1 1 ipn 1 1 — 
1. The inequality (4.1) shows that ijjn is bounded in W^''^ norm. By the Rellich 
theorem [35, Theorem 7.22] we can extract a subsequence, still denoted 
converging to a V'oo G H, weakly in W^'"^ and strongly in L^. Equation (4.6) 
shows that the sequence At/^n converges to zero in L^, and (4.1) with i/j replaced 
with ipn — 1pm shows that ipn is Cauchy in the W^''^ norm. Continuity of A and 
Equation (4.6) imply that ^V'oo = 0, and since A has no kernel on H we obtain 
Voo = 0, which contradicts ||'(/'oo||l2 = 1; and the lemma follows. ■ 

Returning to the proof of Theorem 4.1, define lm{A) to be the image of H 
under A. Then Im(^) is a closed subspace of L^, which can be seen as follows: 

Let tpi be any sequence in H such that the sequence Xi = ^"^i converges in 
to Xoo £ L^- The inequality (4.5) shows that ipi is Cauchy in L^, which together 
with the Carding inequality shows that V'i is Cauchy in the W^'^ norm. As H 
is closed, it follows that there exists ■0oo G H such that ipi converges to V'oo in 
the W^''^ norm, and the equality Xoo = Atp^ follows from continuity of A. 

Let (f) £ L'^ he any element of Im(^)-'-, the orthogonal of Im(^); by 
definition we have 

VV'Gi/ j {(p,AiP)dvY =0 . 

The hypothesis (C2) of elliptic regularity implies that (f) G H, so we can use the 
symmetry of A to conclude 

ytpeH J {A(f),tP)dvY =0 . 

Density of H in implies Acp = 0, thus 

Im(A)-L = Ker(^) . (4.7) 

Define A : H ^ Im(A) by Atp = Aip. By the definition of all the objects 
involved the map A is continuous, surjcctive and injective, hence bijcctivc. 
Let K : Im(^) — > H denote its inverse, then K is continuous by the open 
mapping theorem. Let i be the embedding of W^''^{Y) into L'^{Y); we have 
i{H) C Ker(^)-'- which coincides with Im(A) by (4.7). It follows that for all 
X G Im(yl) we have i o K(x) G Im(A), so that ioK defines a map of Im(74) into 
Im(^), which we will denote by K. Now K is continuous and i compact, which 
implies compactness of K. 

We note that Im(^) is a closed subset of the Hilbert space L^, hence a 
Hilbert space with respect to the induced scalar product. The operator K is 
self-adjoint with respect to this scalar product, which can be seen as follows: 
let ijja = K(/)a, 4'a G Im(^), a = 1, 2, thus tpa & H and Aipa = 4>a- We then have 

J {(j)i,K(j)2) dVY = j (^V'l, V'2) dVY = J {tpl,Atl)2) dVY = J {K(f)i,(f)2) dVY , 

as desired. By the spectral theorem for compact self adjoint operators [63] there 
exists a countable L^-orthonormal basis of Im(A) consisting of eigenfunctions 
of if: 
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with eigenvalues fia accumulating only at 0. Since K is invertible we have 
fXa hence 

A(pa = K<t>a , K = l^a^ ■ 

The required basis of is obtained by completing {(pa} with any L^-orthonormal 
basis of the finite dimensional kernel of A. ■ 



Definition 4.3 A is said to satisfy the spectral condition if A is an operator on 
C°° sections of E over Y which is symmetric with respect to the integration 
pairing with measure dvy and inner product {•,•), and there is a countable 
orthonormal basis {4>a}a&A of L'^T{E) consisting of eigenf unctions, 

Ma = K4>a, a G A , (4.8) 

such that the eigenvalues Aq G M, counted as always with multiplicity, have no 
accumulation point in M. 

Corollary 4.4 Suppose Y is a compact manifold without boundary and A : 

H^{Y) L'^iy) is an elliptic system between sections of the bundles E, F, 
which satisfies the conditions (3.4,3.5) of Theorem 3.7. If A = A^ is formally 
self-adjoint (see (3.7)), then A satisfies the spectral condition, Definition 4-3. 

Proof: Take H = H^{Y). Condition (CO) follows from the coefficient bounds 
(3.4) and the inequality (3.3), and condition (CI) is conclusion (3.22) of The- 
orem 3.7, which also provides condition (C2). Finally, (C3) follows from the 
definition (3.7) of the L^-adjoint A\ since integration by parts is permitted in 
H, and (C4) is standard. The conclusions now follow from Theorem 4.1. ■ 

Corollary 4.5 Suppose Y is a compact manifold without boundary and A : 
H^{Y) — L'^{Y) is an elliptic system between sections of the bundles E, F, 
which satisfies the conditions of Theorem 3.7. There are bases (pa £ ^'^i^), 
ipa G L'^{F), a G a, with real numbers Aq, having no accumulation point in M, 
which satisfy 

Acpa = Xalpa, = ^(pa ■ (4.9) 

The fields (pa, ipa are all H^(Y). 

Proof: This follows directly by applying Corollary 4.4 to the formally self- 
adjoint operator 

which acts between sections of the bundle E ® F. ■ 
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5 Boundary Regularity 



In this section we introduce a broad class of boundary conditions which are 
eUiptic in the sense that the Weak-Strong property (C2) can be estabhshed, at 
least for solutions supported near the boundary. When combined with the inte- 
rior regularity results of §3, this will give the Weak-Strong property for compact 
manifolds with boundary (§6), and for a large class of noncompact manifolds 
with compact boundary (§8). The main result is the boundary regularity The- 
orem 5.11, and the primary ingredient in the arguments is the energy identity 
(5.15) cf. [6] and (2.6). 

We consider operators which may be written abstractly in the form 

L = Lo + B = d^ + A + B, (5.1) 

acting on sections of a (real) vector bundle E over Y x I, where y is a compact 
manifold without boundary^, and for some constant 6 > 0, 

I =[0,6]. 

Let E\y = i*E be the pullback bundle over Y, where i : Y ^ Y x I, 
y (y, 0). We assume that A is an operator on sections of E\y which is 
formally self-adjoint with respect to the pairing defined by integration over Y 
with the measure dvy and the real inner product (■, •) on the fibres of £^|y. The 
operator A and the inner product extend naturally to act on sections of E over 
y X /, and we likewise extend the definition of the integration pairing by using 
the product measure dvydx onYxI. Thus, A is x-independent, but we allow 
B to depend upon x. 

We assume also that A satisfies the spectral condition. Definition 4.3, so 
there is a countable index set A and an orthonormal basis {^^'al^gA -^^(-^ly) 
consisting of eigenfunctions 

A(l)a = >^a<l>a, a G A, (5.2) 

with eigenvalues Aq, G M having no accumulation points in M. 

A formally self-adjoint first order elliptic operator with the coefficient con- 
ditions of Theorem 3.7, will satisfy these conditions, by Corollary 4.4. 

Although we have in mind primarily the case where A, B are first order dif- 
ferential (Dirac-type) operators, the results here will be presented in an abstract 
form, because they could be applied more widely. For example, A = —Ay will 
also satisfy the spectral conditions, so the boundary regularity result Theorem 
5.11 may also be applied to the heat equation. 

Wc fix an eigenvalue cutoff parameter k > 0, which is used to partition the 
index set A into 

{a G A, A„ > k}, (5.3) 

{a G A, Ao < -k}, (5.4) 

{a G A, |Aq,| < k} , (5.5) 

^Although many of the arguments of this section may be extended to allow dY ^ 0, this 
would introduce technical complications which are not relevant to the applications we have in 
mind. 



A+ = 
A" = 
A" = 
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and we set A' = A+ U A . It will be useful also to introduce a scale parameter 

00 = max |Aq,| < 1, (5-6) 

which measures the size of the "small" eigenvalues. For example, we could 
choose K to be the smallest nonzero eigenvalue, k = inf^^^o I^^qI) in which case 
^0 = 0. Choosing k appropriately will lead to estimates for L which are uniform 
under perturbations of A which create or destroy small and zero eigenvalues. 
The eigenfunction expansion u = Yla^A '^a^t^a of -u G L?{Y), where 

Ua ■■= j> {u,(j)a) dvY, 

leads to projection operators Pj^,P-,Pq,P\ defined by 

P±U = ^ UaCpa, (5.7) 
PqU = ^ Ua(i)a , (5.8) 

and P' = 1- Po = P+ + P-. 

For s > 0, the Sobolev-type space H^{Y) is defined as the completion of the 
space of smooth sections C°°{Y), with respect to the norm 

aSA' aSAO 

The space Hl{Y x /) is likewise the completion of C°°{Y x /) with respect to 
the norm 



\q;6A aeA' 



aeA° 



= J j) {\d^u\'^ + \Au\^ + k'^\Pou\'^) dvydx , (5.10) 

where ' = Of course in the typical case where 4 is a first order elliptic 
operator, these norms will be equivalent to the usual Sobolev norms, defined 
using the Fourier transform. Note that the normalization (5.9) ensures that the 
norm is controlled by the Sobolev norm: 

ll^^llz,2(y) < K"1klkj(y), s>0. 

In addition, this formulation leads simply to a useful trace lemma, stated 
in terms of the parameter £, 

e = K5, (5.11) 

which measures the thickness of the boundary layer Y x [0, 6] in units of k~^. 
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Lemma 5.1 The restriction map ry : u i— > vyu = ^(0, •) from C°°{Y) to 
C°°{Y X I), I = [0,(5], extends to a bounded linear map ry '■ Hl(Y x /) — 
i?y^(y) satisfying 

||ryu||^i/2 < ci||n||^i(yx^) , (5.12) 



where c\ = c\{^) = £ "^(1 + \/l + i'^). The m,a,p x ^ ry^xU (where ry^xU = 
u{x,-) is the restriction to Y x {x}), is likewise bounded and continuous in 

1/2 

X from I to Hj {Y). Moreover, ry is surjective: there is an extension map 
ey : Hy^{Y) Hl{Y x I) such that ryeyu = u for all u G hI^'^{Y) and ey 
satisfies ry^sey{u) = and 

Proof: For x e [0, oo) set x{x) = max(0, 1 - x). For any u G C°°{Y x [0,6]) 
with Ua{x), a e A denoting the spectral coefficients, and with x(x) = x(x/(5), 
we find that 

|U„(0)P = -2 / {xUa,:^{xUa))dx 

Jo 

< / ((2x1x1 + vx") Kf + V'^fKf) dx 
Jo 

for any 77 > 0. For a G A' we take r] = a\Xa\ and with an appropriate choice of 
X we find 

\\a\\Ua{0)f < / K'KI" + Mia + 2/£)\Uaf) dx . 

Jo 

Likewise for a G A'^, setting r] = an gives 

Jo 

Choosing a = ^~^(vT+~£2 _ 1) and combining the two estimates gives (5.12) 
for all u G C°°{Y x [0, S]). But this space is dense in Hl{Y x I) by definition, 
and it follows easily that ryu is defined and (5.12) is valid for all u G Hl{Y x I). 
A very similar argument shows that for x G [0,5], 

||ry,.u||J^y2 < ^"'(2 + V4+^)||«||^i(y^j). 

1/2 

To establish continuity of x Ty^xU as a map [0,5] ^ m'%Y), note first 
that for any v e Hl{Y x [xo,xi]), xq < xi, the spectral coefficients lie in 
H^{[xo,xi]) and we may compute: 

\\Va{xi)\'^ - \Vaixo)\'^\ < 2/ \{Va,v'a)\dx 

Jxo 

< Va' [ '(Kl' + ^al^al') dx . 
Jxn 
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Choosing t/q, = |Aq.| for a G A' and rja = k, for a E A° and summing gives 

rY,xiV\\lji/2 - \\rY,xoV\\]ji/2 < I + XI I'^aPl^^aP + ^ K^l^^aN dx 

< lkll//i(yx[2:„,:ri]) • (5-14) 

Given u G Hl{Y x /) and x G [0,6), e G (0, (5 - x)/2), we set 

= u{x + e + x) — u{x + e — x), 

where x G [— e, e] and the ^-dependence of u,v is understood. Applying (5.14) 
with xq = 0, Xi = e, gives v{xo) = 0, v{xi) = u{x + 2e) — n(a;) and 



< 2\\u 



|2 

|2 



Hl(Yx[S,S+2e]) 

= o(l) as e \ . 
This estabhshes continuity from the right, and left continuity follows similarly. 

1 /2 

To see that ry is surjective, we construct an extension map ey : H^' (Y) —>■ 
Hl{Y X /), such that ry o ey = Id. For any u G hI^'^{Y) let {tt^} be an 
approximating Cauchy sequence of smooth fields with spectral coefficients 
and consider the sequence {u'^} defined by 



aGA 



where rja = V^/\\a\ for a G A' and ^/3/K for a G A°, so 

\\i^''\\Hi(Yxi) < [ Y.\^a\''\ui\^{x\x/r]a) + lx'\x/ria))dx 

aeA' 

+ [ Yl '^^K\\x\x/7]o,) + ^x\x/7]o,)') dx . 

Using the bounds x^{x)dx < 1/3, }(''^{x)dx < 1, and noting that 

rd poo 

/ '4)^{x/r})dx < rj '4)^{x)dx 
Jo Jq 

for any ip, we have 

||^i'||^i(yx/) < |Aa||n^|'+ X 

^ aSA' aSAo 

Hence the sequence {u^} is uniformly bounded, and a similar argument shows 
that it is also Cauchy, with limit u = eyu G Hl{YxI). It follows easily that the 
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sequence has boundary values converging to u G hI^^{Y), and u, eyu satisfy 
the bound (5.13). ■ 

The next result relates estimates to boundary conditions, and is the 
key to understanding the natiirc of cllipticity for boundary data. It may be 
considered as a generalization cither of the integration by parts formula for the 
Dirac operator (2.6), or of the estimate underlying the analysis in [6]. We use 
u{0),u{S) to denote the restrictions ryu = uly^^^Qy, ry^su = '"|yx{<5} respec- 
tively. 

Lemma 5.2 Suppose f ^ L?{Y x I) and u e Hl{Y x I) satisfies Lqu = f, then 

IMhUyxI) ^ m-Po)f\\h(y.i) + {i + Olmf\\hiy.i) 
+ 3k'^\\Pou\\12(^YxI) 

+ ^ K{\uam'-K{S)\') 

aeA+ 

+ J2 \>^c.\ {\uam^ - Km^) (5.15) 

aSA- 

Proof: The coefficient functions Ua{x) are measurable and, by Fubini's theo- 
rem, square-integrable over [0,5]. Testing the weak formulation with <j){x,y) = 
x{x)4>a{y) where x ^ C^{{0,6)) shows that Uq, satisfies 



/ 

Jo 



s 

{-UaX' + UaKX - faX) dx = (5.16) 



for all X € C^{{0,S)). Because u G Hl{Y x /), the spectral coefficient Ua{x) 
is differentiable for a.e. x G [0,6], with square-integrable, and (5.16) shows 
that it satisfies the ordinary differential equation u'^{x) + XaUa{x) = fa{x)- 
The trace lemma also shows that the restrictions Ua{0), Ua{5) are well defined. 
Prom the ODE we derive the fundamental identity 

/ \fa\^dx = / \u'^ + XaUal'^ dx 

Jo Jo 

+ -^akaP) dx 

+ Aa(K(^)|2-K(0)|2). (5.17) 

Summing over a G U A~ and noting that the boundary restrictions ■u(O), 
u{6) are in hI^^{Y) by Lemma 5.1 since u G Hl{Y x I) by assumption, we find 



\P'^\\hUYxI) = f (l<l' + |AalVa|')da^ 

° aGA+UA- 

= j> \P'f?dvydx 



+ Yl A«(K(o)l'-KWn ■ (5-18) 

aeA+UA- 
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For a e we use Uq, — fa — AqUq. to estimate 



\\Pou\\hi(YxI) = / ^ iKl'^ + i^'^l'f^al'^) dx 



oeAo 



< / {{1 + s)\U\' + {k' + {1 + s-')4)K\') 

< 3At2||Pon||i.(y,,) + (l + e2)||Po/|li2(yx/) > 



(5.19) 



having chosen e = 6q. Combining (5.18) and (5.19) gives (5.15). 



Remark 5.3 The above proof could be generahzed to allow u ^ LP' and to 
show then that u is in H^^^, but this refinement is unnecessary as we soon will 
show a more general regularity theorem. Working with u G allows us to use 
the boundary terms with impunity — a freedom that is not possible with weak 
solutions at this stage. 

The fundamental estimate (5.15) shows that in order to obtain a useful a 
priori elliptic estimate for a general solution of Lqu = /, it is necessary to 
impose boundary conditions which control P+m(0) (and P-u{5)). Motivated by 
the examples of the spectral and pointwise boundary conditions for the Dirac 
equation (see §2), we introduce a class of boundary conditions which allow 
us to exploit the "good" terms in P_n(0) in (5.15) to provide the required 
control. The effect of the parameterization below is to describe the class of 
admissible boundary data as graphs over the complementary subspace of "good" 
data (1 — P+)ijy^(y). The first justification of this approach is the following 
existence result and its corresponding elliptic estimate (5.24). 

Lemma 5.4 Let P = -P^^+uA spectral projection determined by A'^ and 

some subset A C A° o/ the set of small eigenvalues. Let a G PLlI^'^{Y) and 
f G L2(y X I) be given, and suppose K : {1 - P)hV'^{Y) PH^^iY) is 
a continuous linear operator, so there is a constant k > such that for all 



Then there exists a solution u G Hi (Y x /) to the boundary value problem 



w G hI^\y) 



\\K{l-P)w\\„^/2 <k\\{l-P)w\\„^/2. 



(5.20) 



LqU 

Pu{0) 

(i-PHS) 



f 

a + K{l-P)u{0) 
0. 



(5.21) 
(5.22) 
(5.23) 



Moreover, the solution u satisfies the estimate 




2 

Hl{YxI) 



) 



(5.24) 



where the constant C4 depends on i,6o and k. 
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Proof: The solution to an ordinary differential equation u'{x) + Xu = f may 
be written in either of the two forms 



u{x) = < 



J\s (5.25) 



J X 



Consider first the spectral coefficients Ua{x) for a G A~ U A', where A' = A''\A. 
The boundary condition (5.23) is achieved for Ua by letting Ua{S) = 0, so we 
define 

u^(x) = - e^«(^-^Va(s) ds, aeA~U A', (5.26) 

J X 

where faix) is the spectral coefficient of /. Note that fa £ L'^i^), so the integral 
in (5.26) is well-defined. The identity (5.17) and «q(5) = shows that 

[\\u'af + Xl\ua\'')dx = Xa\uam''+[^\fa\''dx, Va € A" U A' . (5.27) 
Jo Jo 

It follows that for all a G A~, 

|Aa||^^a(0)P < r Ifal^dx - + Xl\Ua\'') dx. (5.28) 

Jo Jo 

To control the small eigenvalues a G A' we use an elementary lemma, the 
proof of which is an exercise: 



ai: 



e^('-^)f{s)ds] dx<{ •^o-' (5.29) 



Lemma 5.5 For any f e i ([0, S]) and XeR, 

\5^e^^^ f{x)dx A>0 
^d'^J^f{x)dx A<0. 

Prom (5.26) and Lemma 5.5 it follows that for a G A', 

s s 

[ K^al^dx < i^V^^o / \fafdx. 

Jo Jo 
Using u'^ = fa — ^a^a ^ (5.19) we obtain 

/ i\u'a\^ + K^K\^)dx < {1 + C2) [ Ifal^dx, (5.30) 

Jo Jo 



where 

C2 = C2(A0o) = ^o + i^V^^°- (5-31) 
Combining (5.27) and (5.30) shows that 



{x,y) := ^ Ua{x)(j)a{y) 
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is a sum converging in Hl(Y x /), and ^ satisfies 

h^'^WHliVxI) ^ \\i^-P)f\\h(Y.I)+C2\\PAJ\\h(Y.I) 

- Yl l^«||^^a(0)|' , (5.32) 

aeA- 

where 1 — P = P_ + Pj^,. For a G A' such that Aq, ^ 0, using (5.26) we find 
that 

2 



< ^(e''"'-l) Tl/. 

^^a Jo 



5 

a I 



while if Aq = then |nQ,(0)p < S \ fa\'^dx. Defining 

C3 = £e'^^^°, (5.33) 

it follows that 

K|Ma(0)|^ < C3 / Ifal'^dx , V « G A'. 

Jo 

Thus, combining with (5.32) we have 

Il^^-^ll|i(yx7) + 11^^-^0)115^,/, 

< 11(1 - P)f\\hiYxi) + (C2 + C3)||PA'/lli2(yx7)- (5-34) 
For a G A+ U A we use (5.25) to define Ua by 

Ua{x) = e-^"^(aa + Kau(-)(0)) + re^"(^-^)/„(s)ds, (5.35) 

Jo 

where (Jq, KceU^''\0) denote the (pa coefficients of a and ^-^("^(0) respectively. 
Note in particular that (5.34) shows that u(-)(0) G hI^^{Y), so Ku^-\0) G 
Hy^{Y) by the hypothesis (5.20), hence the coefficients Kau'^''\0) are well 
defined. 

For a G A+ we estimate using (5.17) and (5.22): 

s s 
f (Kl' + A^)dx < A,K + i^,u(-)(0)p+ / \fafdx 
Jo Jo 

1-5 

< 2\a{\aa\^ + \Kau(-\0)\^)+ \fa\^dx . 

Jo 

For a G A we estimate directly from (5.35): 

\ua\'^dx < SK^aal"^ + \KaU^-\0)\'^) f\-^^"^dx 
Jo Jo 

+ 3^2 e^'^(*-^)/a(s)ds j dx 

< 3c3Ki\aaf + \KaU^-\0)f) + f^Cs / \fa\^dx , (5.37) 

Jo 
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where Lemma 5.5 has been used to control the final term. Using u'^ = fa — ^aUa 
to estimate \u'^\'^ < {l + e)\fa\'^ + {1 + e-'^)Xl\ua\'^ with e = 9^, (5.37) gives for 
a e A, 

/ (|<P + k2|„^|2)^^ < 9c3K(|aa|2 + |Ke,uH(0)|2) 

+ (1 + 3C2) / \fa?dx. (5.38) 
Jo 

Combining (5.36) and (5.38) we have (setting ■u(+) = X^A+yX Ua4>a) 



\u 



+ 2{\\P+a\\ly, + \\P+Ku^-\0)\\ly.) 

+ 9csi\\P^a\\l,j, + \\P^Ku^-\0)\\ly,) , (5.39) 



where P+ = -Pa+, P = P^ + P^. Since we have already shown that KaU^ \0) G 

1/2 

(y), all terms on the right hand side of (5.39) are bounded, which shows 
that u^^^ is well-defined in Hl{Y x /). 

With u = u^'^^ + u^^\ we add an appropriate multiple of (5.34) to (5.39) to 
control the bad terms in Ku^~\0) with the good term n^~)(0) of (5.34). This 
gives the elliptic estimate (5.24): 



\u 



\HliYxi) ^ ll^i^^^ll|i(Yx7) + max(l,2A;^9c3A;2)||n( 
< cl(||/||i.(^,,) + ||a||^,/,), 



where C4 = C4{£,9o,k) as required. The definitions (5.26), (5.35) ensure w is a 
solution satisfying the boundary conditions (5.22), (5.23). ■ 

Explicitly, we may take C4 = 3c2 + (fc^ + 1)(2 + 9C3) in general, and C4 = 
2max(l,A;2) if = 0. 

The next result is the key to handling operators with coefficients depending 
on X. Recall that the operator norm ||-B||op of a linear map B : Xi X2 
between Banach spaces is the smallest constant such that 

\\Bu\\x, < \\B\\op\\u\\x^, V« G Xi . (5.40) 

Lemma 5.6 Suppose Lq, A, f, a, K are as in Lemma 5.4, and suppose B : Hl(Yx 
I) — > L^(y x /) is a linear map satisfying 

C4 \\B\\op < I. (5.41) 
Then there exists u G Hl{Y x /) satisfying 

(Lo + B)u = f 
and the boundary conditions (5. 22), (5. 23), such that 

W^'WhiYxi) ^ 1-c1b\L ("^'I^^(>-x7) + lkll^i/2). (5.42) 
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Proof: Let u^^'' e Hl{Y x I) be any function satisfying the boundary condi- 
tions (5. 22), (5. 23); the trace lemma 5.1 ensures the existence of a suitable u^^\ 
Construct a sequence {u^''^} C Hl{Y x /) by solving the problems 

Lo^W = f-Bu^^-^^ (5.43) 
Pu(^)(0) = (j + K{l-P)u^^\{)) (5.44) 
{l-P)u^^^\5) = 0, n = l,2,.... (5.45) 

Lemma 5.4 ensures this problem has a solution for every n > 1, and the differ- 
ence w^^'^ = M^*^) — u^^~^^ satisfies 

Pw;('=)(0) = K{1- P)w^''\Q) 
{l-P)w^''\6) = 0. 

The elliptic estimate (5.24) gives 

\\w^^^\\Hl[Yy<I) < C4\\Bw^''-^^\l2^yxI)- 

If 1 1 -B 1 1 op < l/c4 then the iteration is a contraction and thus the sequence u^''^ 
is Cauchy, converging to u = lim„_^oo strongly in Hl{Y x I). Taking the 
limit of (5.43) shows that {Lq + B)u = f, and boundedness of the trace operator 
ry shows that u satisfies the boundary conditions (5.22)-(5.23). The elliptic 
estimates (5.24) satisfied by u^''^ are preserved in the limit, so u satisfies 

\\u\\hI{YxI) < CiiWf - Bu\\l'2{YxI) + Ikllj^i/a), 

from which (5.42) follows easily. ■ 

Observe that the proof of Lemma 5.6 relies on just two properties of the 
operator Lq; namely, the solvability of the problem (5.43) with boundary con- 
ditions (5.44), (5.45), and the elliptic estimate (5.24), which provides the size 
bound (5.41) for the perturbation B. This suggest that it should be possible to 
extend this existence result to more general operators L = Lq + B, for which a 
strictly coercive estimate such as (5.24) can be established. 

Consider, for example, the case where E has a complex structure J : E ^ E, 
= —1, and A is a normal operator ([^,y4*] = 0), so A = Aq + JAi where 
Aq,Ai are self-adjoint and commuting, Aq satisfies the spectral condition, and 
both commute with J. Then A admits an eigenfunction basis 

= (Aa + IJ-aJ)(l>a, f^a € I^, Va G A, 

and the results of this section extend with only minor modifications, provided 
the eigenvalues Xa, fJ'a satisfy the sectorial condition 

SUp|/Xa|/(l + |Aa|) < OO. 

Before stating the main uniqueness theorem, a definition of weak solution 
with boundary conditions is required. Note that although the definition is 
consistent with just boundary data, the regularity theorem 5.11 will require 
data in h'J'. 



28 



Definition 5.7 Suppose L = dx + A -\- B where A satisfies the spectral con- 
ditions (Definition 4.3) and B : Hl(Y x I) ^ L^(Y x I) is a bounded linear 
operator for which there exists an L"^ -adjoint B'^ : Hl{Y x I) —>■ L'^{Y x I) such 
that: 

/ {Bu, v) dvydx = / {u, B^v) dvydx , y u,v e hI{Y x I). (5.46) 

JyxI JyxI 

Suppose further that P = P-\. + P^ (as in Lemma 5.4), that K : (1 — P)L^(y)— > 

PL^(y) is a hounded linear map with L'^ -adjoint K'f : PL'^{Y) {l-P)L'^{Y), 
and let a G PL'^{Y), f e L^{Y x I) be given. A weak solution of the boundary 
value problem 

Lu = f (5.47) 

Pu{0) = a + K(l - P)u(O) (5.48) 

{l-P)u{S) = (5.49) 

is a field u G L^(y x /) satisfying (with = —dx -\- A-\- B^) 

[ {u, L^(f)) dvydx = I (/, (j)) dvydx + / {a, 4>{0)) dvy , (5.50) 
Jyxi Jyxi Jy 

for all 4> G HI{Y x I) satisfying the adjoint boundary conditions 

(1 - P + Ktp)0(O) = (5.51) 
P(j){6) = 0. (5.52) 

1/2 

The boundary values 0(0), ^(1) both lie in H^, (Y) by the trace lemma, 
so the adjoint boundary conditions are well-defined on the space of test fields. 
Since C°° fields are dense in Hl{Y x I) and in Hy^{Y), to verify the weak 
equation (5.50) it suffices to test just with C°° fields (j); however the uniqueness 
argument of Lemma 5.10 requires the use of an Hi test field. 

The structure of the adjoint boundary condition (5.51) is explained by the 
next lemma, which is applied with v = u{0) — a and H = L'^{Y). 

Lemma 5.8 If H is a Hilbert space, P : H ^ H is an orthogonal projection 
and K : ker P — ^ range P is bounded, and if v E H satisfies 

{v,(t))H = V G ker(l - P + i^tp) (5.53) 

(where is the adjoint of K in H), then v G ker(P — K{1 — P)). 

Proof: Since ker P = range (1 — P) _L range P, it follows that P is self-adjoint 
and there is an orthogonal decomposition H = {1 — P)H © PH. Setting 
(Pi = (1 - P)(P, (j)2 = P(j), the condition = </)2 G ker(l - P + K^P) 

is equivalent to (pi = —K'^(p2, which exhibits ker(l — P + K'^P) as a graph over 
PH. Similarly decomposing v = vi -\- V2, the condition {v, (p) = is equivalent 
to {v2 — Kvi,(p2) = 0. Since this holds for all (p2 G PH, it follows that V2 = Kvi, 
or equivalently, v G ker(P — K{1 — P)). ■ 
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In other words, H = ker(P — K{1 — P)) © kcr(l — P + i^^P) is an orthogonal 
sphtting of H, where P — K{1 — P), 1 — P + K'^P are projections, which are 
not orthogonal in general. 

Using Lemma 5.8 we next show that an weak solution of the bound- 
ary value problem (5.47)-(5.49), in fact satisfies the equation (5.47) and the 
boundary conditions (5.48,5.49) in the strong sense: 

Lemma 5.9 If u ^ Hl{Y x I) is a weak solution of (5.47) with the boundary 
conditions (5.48), (5.49), then u satisfies the equation Lu = f in the sense of 
strong (H^) derivatives, and the restrictions u{0) = ry(«), u{5) = ry^su satisfy 
the boundary conditions (5. 48), (5. 49) in L'^iY). Conversely, if u E Hl{Y x I) 
is a strong solution of (5.47,5.48,5.49), then u is also a weak solution. 

Proof: Integration by parts gives 



for any u,(j) e Hl(V x I). Testing u with arbitrary ^ G C^{Y x 7) shows that a 
Hi weak solution satisfies Lu = f in the sense of strong derivatives. Comparing 
this formula with (5.50) shows also that 



for all 0(0), (t>{5) satisfying the adjoint boundary conditions (5.51), (5.52). Since 
Hy^iY) is dense in ker(l - P + K'fP) C L'^{Y), Lemma 5.8 may be applied 
with V = u{0) — (7 to show the boundary condition (5.48) holds in L'^{Y), and 
(5.49) follows similarly. 

To show the converse, integration by parts again gives 



and the final two terms vanish by the adjoint boundary conditions (5.51,5.52). 
■ 

By solving an adjoint problem, we now show that weak solutions of (5.47)- 
(5.49) are unique. 

Lemma 5.10 Letu G L'^{YxI) be a weak solution of the boundary value problem 
(5.47)-(5.49), with a G L'^{Y) and f £ L'^{Y x /). Suppose that the operator 
L = dx + A + B satisfies the conditions of Lemma 5.6 and Definition 5.7, and 
the L^{Y X L) -adjoint Pt ; hI{Y x I) ^ L'^{Y x I) satisfies 







{a, 0(0)) + (u(0), (1 - P + K^P)0(O)) - («(<5), P0(5)) , 



C4||-B'''||op < 1- 



(5.54) 
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Suppose also that the boundary operators K, of Definition 5. 7 satisfy 

\\K{l-P)w\\^,,^ < k\\{l-P)w\\^,/2 , (5.55) 
\\K^Pw\\^i/2 < k\\Pw\\^i/2. (5.56) 

1/2 

for some constant k > and all w & Hj (Y). Then u is unique. 

Proof: It will suffice to show that any weak solution u of (5.47)-(5.49) with a = 
0, / = 0, must vanish. Consider the adjoint problem L^<p = u with boundary 
conditions (5. 51), (5. 52); writing L'^cj) = u as {dx — A — B'^)(f) = —u, we see 
that L^, satisfy the conditions required by Lemma 5.6, since interchanging 
A ^ —A means replacing P by 1 — P, and perhaps changing a finite number of 
eigenfunctions in A (without modifying Aq). The elliptic estimate (5.42) does 
not depend on A. Thus, by Lemma 5.6 there exists a solution (p G Hliy x /) 
of this boundary value problem. By construction, ^ satisfies the boundary 
conditions required of test functions in (5.50), so testing u in (5.50) with 
gives 

/ = / {u,L^4>) = 

JyxI JyxI 

and thus u = 0. ■ 



It is easy to check that (5.56) is equivalent to requiring that K : (1 — 
P)H~'^^^{Y) PH:;^''^{Y) is bounded, with constant k. 

We now obtain the main result on boundary regularity of weak solutions. 
Note that although the definition of weak solution assumes boundary data a G 
L^(y) only, and uniqueness of weak solutions holds also in this generality, this 
condition is incompatible with regularity u G .ff^ (K x /) , which would imply (by 
simple restriction) that a G h\^'^{Y). However, some results for boundary 
conditions on domains with uniformly Lipschitz boundary are known [1,2], so 
it is plausible that the results here could be extended. 

Theorem 5.11 Suppose u E L'^{Y x I) is a weak solution of the boundary 

value problem (5.47)-(5.49) with operator L = dx + A + Bq + Bi, where A 
satisfies the spectral conditions (Definition 4-3), Bq satisfies the size condition 
(5.41) with L'^-adjoint B^ satisfying (5.54), and Bi : L'^{Y x I) ^ L'^{Y x I) 
is bounded. Further suppose the boundary operators K,K^ satisfy (5.55,5.56), 
and a G PhV^{Y). Then u G Hl{Y x I) (so u is a strong solution) and u 
satisfies the a priori estimate 

Il'"llifi(yx7) < 1 _ ^ ifp II (ll/IU2(rxJ) + + ||-Bl||L2-»L2||^*llL2(yx7))- 

i C4||J3o||op * 

(5.57) 

Proof: Since ||i?i?x||^2(yx/) < ||-Bi||op||?i||L2(yx/)) ^ satisfies Lu := {dx + A + 
Bo)u = f where f = f — Biu e L'^{Y x I). Lemma 5.6 constructs a solution 
u G X I) of Lu = f satisfying the same boundary conditions, and it 

follows that u is also a weak solution. By the Uniqueness Lemma 5.10 we have 
u = u and thus u G Hl{Y x I), as required. The estimate (5.42) of Lemma 5.6 
(with / replaced by /) leads directly to (5.57). ■ 
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6 Boundary regularity for first order systems 



In this section we determine conditions under which the boundary regularity 
results of §5 apply to a first order equation of Dirac type (see (6.7), (6.9)) at 
the boundary, with suitable boundary operator, to show regularity of an 
weak solution. 

We assume M is a smooth manifold with compact boundary Y, E M 
and F ^ M are real vector bundles over M with scalar products, and C is 
a first order elliptic operator on sections of E to sections of F, which in local 
coordinates , j = 1, . . . ,n takes the form 

Cu = d'djU + bu , (6.1) 

where , b are homomorphisms of £^ to F as before. Note that we do not 
assume F to be connected. 

To apply the preceding results, the coefficients must satisfy the interior 
regularity conditions (3.4), and the boundary restrictions must be defined and 
satisfy the corresponding conditions in dimension n — 1, 

a^\y e W^i'("-i)*(y)nCO, j = l,...,n, 

b\y e L(--'r{Y). ^^-^^ 

Conditions (3.4,3.5) and (6.2) will be assumed throughout this section. 
Remark 6.1 The conditions 

G <f(M)nCO(M), 

b e W(-'^\M), ^'-'^ 

where 

s = n/2 for n > 4 , 

s > 3/2 for n = 3 , (6.4) 
s = 3/2 for n = 2, 

imply the interior (3.4) and boundary (6.2) coefficient regularity conditions, 
through the Sobolev embedding and trace theorems [60]. The C° condition in 
(6.3) is superfluous for n = 2, 3. 

Let x = be a boundary coordinate, defining a tubular neighbourhood 
y X [0, 1] C M of F with local coordinates {y\x) € F x [0,1] (where we 
identify Y with Y x {0}). Let dvM be a volume measure on M, and define 
dvy = {—l)"'~^dx JdvM\Y on y x {0}. The local coordinate integration factor 
7 is defined in y x [0, 1] by d,VM = 'y dy dx, where dy dx is coordinate Lebesgue 
measure. We assume the local coordinate condition^ 

7G (w^'"*nc°)(y X [0,1]) . (6.5) 

In order to directly apply the results of the previous section, we assume that 
7 = 7(0) is independent of a; in y x [0,1], so dvM = dvydx. This involves 

^As in §3, this means that we can cover K by a finite number of coordinate charts ffa so 
that 7 has the stated regulaxity in the local coordinates on 0'a x [0, 1]. 
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no loss of generality, as the x-dependence of 7 in the integral form (6.15) of 
the equation can be absorbed through a rescaling of the coefficients a*, b. Since 
the restrictions to y x {0} are unchanged, this does not affect the boundary 
operator A. 

To minimize confusion with the outer unit normal n = —dx we set 

u := *a" , (6.6) 

and for simplicity we assume 

V = u-^ . (6.7) 

This is satisfied by the Dirac operator. Equation (6.7) can also be achieved in 
several other situations of interest by pre-multiplying £ by a suitable homo- 
morphism, or by making a frame change in F (this can be done e.g. when F is 
trivial). Note that such a pre- multiplication does not affect the values of a*, b, 
where we define 

a' := z^a*|y, i = 1, . . . ,n — 1, b := vbly ■ (6-8) 

By extending independent of x we regard a*, 6 as defined on Y x [0,1]. We 
assume the important boundary symmetry condition 

a' = - *a\ i = l,...,n - 1 , (6.9) 

where the transpose is taken with respect to the inner product on E. We 
wrap up these conditions into a definition. 

Definition 6.2 A first order system (6.1) is of boundary Dirac type if the 
coefficients satisfy the conditions (6.7,6.9) in some neighbourhood of the 
boundary. 

As discussed in §10, this class includes the examples of §2 and §11. These 
conditions will be assumed henceforth. 

Using (6.9) we define the boundary operator 

n-1 

An := a^djU + apu + bpu , (6.10) 

1=1 

where 

n-1 

ao = lJ2{dia' + a'dilogj) , (6.11) 

i=l 

and bo = 60 |y for some symmetric endomorphism bo of E. We require that 
bo G L"*(y X /)) and 60 G L("-i)*(y), compare (3.4), (6.2). Then A is formally 
self-adjoint {A^ = A) on the bundle E\y over Y, with respect to the measure 
dvy. 

Note that the choice of zero-order term bo gives some freedom in the defini- 
tion of the boundary operator A, which is thus not uniquely determined by jC. 
Near the boundary, Cu = ^f^dx -\- A + B), which may be expressed as 

Lu:={dx + A + B)u = f , (6.12) 
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where Lu = vCu, f = vf, and B denotes the difference 

Ti-l 



Bu = ^{i^a' - a')diU + {uh -clq- bo)u . (6.13) 

i=l 

By Corollary 4.4, under the above conditions, A will satisfy the spectral 
condition 4.3. Denote the eigenvalue index set by A and fix a cutoff k, as in §5. 
Similarly let A+ = {a G A : Aq, > «}, fix some subset A C {a G A : |Aa| < k} 
and let P = Pj^+ + be the associated spectral projection operator. Note 
that Theorem 3.7 (see (3.22)) shows that the norms defined using A will be 
equivalent to the corresponding norms defined on Y and y x I, at least for 
s = 0, 1; it seems likely that this will hold, by interpolation, for all s G [0, 1]. 
If the coefficient regularity allows an H'''^^ elliptic estimate (c/. Theorem 3.8) 
then this should extend to s G [0, A; + 1] . 

The definition of weak solution on a manifold with boundary which we are 
about to give is slightly simpler than the tubular neighbourhood Definition 5.7 
of §5, since the conditions (5.49), (5.52) may be imposed by localising with a 
boundary cutoff function; see the proof of Theorem 6.4. Let H^{M) denote the 
dense subspace of H^{M) consisting of functions of compact support, where 
we recall that because M is a manifold with boundary Y = dM, H^{M) in- 
cludes functions which are non-zero on y. As in §5, the boundary condition 
is expressed using a positive spectrum projection P : L'^{Y) L'^O^) and a 
bounded linear map K : {1 - P)L'^{Y) -> PL'^{Y) and its L'^{Y) adjoint . 

Definition 6.3 Let f G Lf^^iM) ^ PL'^{Y) be given. A section u G 

L^Qj,(M) is a weak solution of Cu = f with boundary condition 

Puo = a + K{1 - P)uQ , (6.14) 

if 

/ {u,£^(t>) dvM = {f,(j))dvM+f{(^,i^4>o)dvY, (6.15) 

JM JM JY 

for all 4> G H^{M) satisfying the boundary condition 

{l-P + K^P){i,cf>o) = , (6.16) 
where is the adjoint given by (3.7). 

Note we are using the notation UQ,(j)Q, etc., to denote the restriction (trace) on 
the boundary Y. The additional term in (6.15,6.16) (cf. (5.50,5.51)) arises 
from the relation C = ^^{dx + A + B) between C and the boundary form 
dx + A + B used in §5. 

The boundary condition (6.14) restricts uq = u\y to lie in the affine subspace 
of L'^iY) given by the graph of x i-^ a+Kx over the negative spectrum subspace 
X G (1 — P)L^(y). It will be useful to re-express (6.14) as 

/Cuo = cr , (6-17) 
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where we have introduced the operator /C on L^(y), 

K:=P- K{1 - P) , (6.18) 
and hkewise to re-express the "adjoint" boundary condition (6.16) as 

/CVo = , (6.19) 

where we define 

■=v-^{l-P + K^P)v . (6.20) 

The next result generahzes the interior weak-strong Theorem 3.7 to bound- 
ary value problems. 

Theorem 6.4 Suppose C and A satisfy the conditions (3.4,3.5,6.2, 6.5,6.7,6.9), 
and suppose a G PhI^'^^Y), f e Ll^iM). PuHher suppose K : {l-P)L'^{Y) 
PL?(Y) is bounded linear and satisfies (5.55), with adjoint satisfying 
(5.56). Assume u £ L^^d^) 

is a weak solution of £u = f with the boundary 
condition ICuq = a (6.17). Then u G H^^^{M) and u is a strong solution. 
Moreover, there are constants S G (0,1], C5, depending only on K,k and a^,b, 
and intervals I' = [0,6/2], I = [0,6] such that 

\\u\\m(Yxr) < C5(||/||L2(yx/) + \W\\hI^\y) + \MlHYxI)) ■ (6-21) 

Proof: Theorem 3.7 ensures u G H^^^{M), where M is the interior of M, 
so it suffices to consider u compactly supported in K x [0,6), for any choice 
of 5 G (0,1). In particular, because u then vanishes near Y x {6}, it follows 
from (6.14) that ti is a weak solution with boundary conditions, in the sense of 
Definition 5.7. 

It will suffice to show, for a sufficiently small choice of 6 > 0, that we may 
decompose B = Bq + Bi into pieces satisfying the size conditions of Theorem 
5.11. Write B = 13' di + /? where 

f5\y, x) = (a-(y, x))-^a\y, x) - {a^{y, 0))-\\y, 0), i = 1, . . . , n - 1, 

so G VF^'™* n C° and /3*(y, 0) = 0. Since the constant C4 of Lemma 5.4 
depends only on < !> ^ = ^^'6 and the constant k of (5.55,5.56), it is bounded 
uniformly m. 6 < 1. Consequently for any e > there is 5o > such that 
caW^\\l^{Yx[q,s]) < e for all 6 < 6^. 

Likewise, since 7 G W^'"* (6.5), we have /3 G L"* and there is a decom- 
position /3 = /?o + /9i with C4 C5 ||/?o|Il"* (yx[o 1]) — ^' where Cs is the Sobolev 
constant on y x [0, 1], and /3i G Then Bq = f3'^di + f3o satisfies (5.41), as 
does Bq (possibly after decreasing 6), and Bi = (3i is bounded on L^. Theorem 
5.11 now applies and shows u e H^{Y x [0,6]), since H^{Y x I) = Hl{Y x I) 
as remarked above. The elliptic estimate (6.21) follows by applying (5.57) to 
u = xu, where x = x(^) is a cutoff function, x(a;) = 1 for < x < 6/2, xi^) = 
for X > j6. ■ 
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Corollary 6.5 Suppose A, K, K, satisfy the conditions of Theorem 6.4- Then 
for all u G H^^^{M) we have the boundary estimate 

Il'"lli?i(yx7') < C5{\\£-u\\l2^yxI) + \\f^uo\\j^i/2^y^ + ||'u||L2(yx7)) • (6-22) 

Proof: Hue Hl^{M) then f := Cu e Ll^{M), a := /Cuq G ^^^(y), and 
M is a strong solution. Since integration by parts may be applied to show u 
satisfies the weak equation (6.15), Theorem 6.4 applies and gives (6.22). ■ 

A bootstrap argument, slightly more complicated than that used for the 
interior bounds Theorem 3.8, leads to higher, H^~^^, regularity. Rather than 
stating complicated conditions for general k, we describe the details only for 
the case k = 1 {u ^ H^). The coefficient regularity conditions are most likely 
not optimal. 

Theorem 6.6 In the setting of Theorem 6.4 let u G H^{Y x /) be the solution 
and suppose the following additional regularity conditions are satisfied, 

,b eW^'^{Y X I) , (6.23) 

7 G W^2'"*(y X [0,1]) , (6.24) 

feH^YxI), aeH^/^Y), (6.25) 

[A, K]{1 - P) : (1 - P)hI^^{Y) PHI^^{Y) is bounded. (6.26) 

Then there exists 6" < S/2 such that u e H'^{Y x I"), where I" = [0,6"], and 
there is a constant cq depending on the coefficient bounds (6.23)-(6.26), such 
that 

Il^*lli72(yx/") < CQ{\\f\\m{Yxi) + Il^llij3/2(y) + ll^llL2(yx7)) ■ (6-27) 

Remark 6.7 For the APS and chiral boundary conditions (2.13), (2.18), A 
commutes with K and thus (6.26) is trivially satisfied. 

Proof: The idea is to show that Au satisfies a similar boundary value problem. 
For convenience, let Hq{Y xl), I = [0, (5], denote the completion of the C°° 
functions of compact support in y x [0, 8), where 5 is the constant of Theorem 
6.4. In particular, functions in Hq{Y x I) have vanishing trace on y x {5}. For 
any ti, G C'^{Y x [0, (^)), we have the identity 

/ {Av^Dip) dxdvy = / [{[L,A\v,'^) + {Lv,A'il))\dxdvY 
JyxI JyxI ^ ' 

+ (f {vo,Ai;o)dvY . (6.28) 



Since A is formally self-adjoint with respect to dvY, this formula follows by 

direct calculation. The terms with L'l''0, Lv are well-defined for v^ip ^ Hq{Y xl). 
1/2 

Since fo,^o £ Hj (Y) by Lemma 5.1, the boundary integral extends also by 
writing it as 

j) {vo, ^Vo) dvY = <j> {Jvo, J'-^Aipo) dvY ■ (6.29) 
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Here J = (1 + |^|)^/^ is defined as Ju = X^asA'^all + \^Q\Y^'^i>a, with u = 
J2aeA''^oi'^oi', is defined similarly. This shows that 

j>^{vo,A'4)Q)dvY < c\\vQ\\jji/2^Y)\\'^o\\jji/2^Yy 

1 /2 

for all vojipo G H^' {Y), where the constant c is determined by A. 

More care is required to control the commutator [-L,^], which takes the 
form 

[L, A]v = a'^dfjV + a\diV + a2V 

where 

a\ = va-^ djO,^ — a-' dj{i'a'') + [va^ ,ao], 
02 = ua^diQQ — a^di{vh) . 

Observe that a'-'ly = 0. The coefficient conditions (6.23) and v & ensure 
that a\diV is in and hence may be combined with the source term Af. 
Likewise (6.23) ensures a2V is bounded in L^. 

Let Ao G R\A, so ^ — Aq has trivial kernel and satisfies elliptic estimates 

ll^^lliys+i < c||(^ — Ao)f||_f/s for s = 0, 1 at least. Since A is self-adjoint and 
elliptic, the cokcrncl of ^ — Aq is also trivial, so ^ — Aq : H^~^^ — >■ is 
invertible. Now decompose a^^df^v = B2AV + B^v where 

S2=a'^'4(A-Ao)-^ 
is bounded from L^, and 

Bs = -Xoa'^dfj{A-Xo)-^ 

is also bounded from L? . Note that by perhaps decreasing 5 we may 

ensure that B2 and B2 satisfy a smallness condition similar to (5.41). 

Direct calculation (noting that [A, P\ = 0) establishes the boundary formula 

j) {{1 + K){1- P)vo,A2Po)dvY = ^ {Avo,{l-P + K^P)ipo)dvY 

+ - P)vo,i^o) dvY i^.SO) 

for all vojipo G hI^^{Y), since K satisfies (6.26) and (5.55) by assumption. 

Now u G Hq{Y X I) satisfies Lu = f and uq = cr + (1 + K)(l — P)uo. 
Substituting u for v in (6.28) and using these relations, shows that u satisfies 

/ {Au, {L — B2)^ip) dx dvy = / {Af + a\diU + a2U, ijj) dx dvy 
JyxI JyxI 

+ £{Aa + [A, K]{1 - P)uQ, 1P0) dvy 
+ j {AuQ, (1 - P + K^P)il)o) dvY{£>M) 
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In particular, if Vo e ker(l -P + K^P) fl hI^^{Y) then w = Au e L'^{Y x I) 
satisfies 

/ {w, {L — 82)^ ip) dx dvY = / (/i, ■0) dx dvY 
JyxI JyxI 

+ £{ai,2Po)dvY, (6.32) 

for all £ H^(Y X I) such that ip £ kei{l - P + K^^ P) , where 

fl = Af + a\diU + a2U £ , 

(71 = Aa + [A, K](1-P)uq£ hI'"^ {Y) . 

In other words, G L^(y x /) is a weak solution of the problem 

{L-B2)w = fl, 

Pwo = ai + K{1 - P)wo . 

By shrinking the boundary layer we may assume ||-B2||ifi^L2 and ||5|||^i_>^2 
are sufficiently small that the conditions of Theorem 6.4 are met, so w = Au G 
H^{Y X [0,(5]). The equation now gives dxU = f — Au — Bu G and thus 
u£H^{Y x[Q,6]). ■ 



7 Fredholm properties on compact manifolds 

The interior and boundary estimates of § 6 lead to solvability (Fredholm) results, 
by standard arguments. The main interest lies in identifying the cokernel, and 
we give a simple necessary and sufficient condition for solvability, in Theorem 
7.3. This section treats only compact manifolds, leaving the more difficult 
case of non-compact manifolds to the following section. Because more detailed 
descriptions are given in § 8, some of the arguments are only briefly summarised 
here. 

Throughout this section we assume the coefficients a-^ , j = 1, . . . , n and h of 
C satisfy the conditions of § 6, namely (3.4), (3.5), (6.2), (6.5), (6.7), (6.9); 
is given by (3.7), and the boundary operators K. K'^ satisfy (5.55,5.56), where 
P = Ply + Pf^ is a positive spectrum projection of A (6.10), and 1C ,lO are 
defined by (6.18,6.20). 

Recall the Sobolev space H^{M) of sections of E over M is defined by the 
norm (3.17) 

\MhHm)= [ i\Vuf + \uf)dvM, (7.1) 

where lengths are measured using the metric ( , ) on E and a fixed smooth 
background metric g on TM, and the connection V satisfies (3.18,3.19). Note 
again that V need not be compatible with the metric on E, and the space 
H^{M) is independent of the choice of V. 

The following basic elliptic estimate extends (3.22) of Theorem 3.7 to man- 
ifolds with boundary, using the boundary neighbourhood estimate (6.21) of 
Theorem 6.4. 
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Proposition 7.1 There is a constant C > depending on ,b,T and K such 
that for allue H^{M), 

\\u\\m{M) < C{\\£.u\\l2(^m) + \\J<^uo\\^i/2^y^ + \\u\\l2(^m)) ■ (7.2) 

Proof: The argument used in Theorem 3.7 to prove the interior estimate (3.22) 
may be appHed using Theorem 6.4, estimate (6.21), to estimate ||^ta|| over 
boundary neighbourhoods Ua- The remaining details are unchanged. ■ 

Theorem 7.2 The linear operator 

(£, /C) : (M) ^ (M) x PhI^^ {Y) (7.3) 

is semi-Predholm (i.e. has finite dimensional kernel and closed range). 

Proof: Suppose {w,fc}i° is a sequence in ker(£,/C), normalised by HufcHji^ii-j;^) = 
1. To show the kernel is finite dimensional, it suffices to show there is a subse- 
quence converging in (M) . By Rellich's lemma there is a subsequence (which 
we also denote Uk) which converges strongly in L'^(M), to u £ L'^{M) say. The 
elliptic estimate applied to the differences Uj — shows that the sequence is 
Cauchy in H^{M) and thus converges strongly to u & H^{M). Since (7.3) is 
bounded, it follows that u G ker(>C, /C), so the unit ball in the kernel is compact 
and hence the kernel is finite dimensional. 

To show the range is closed, let H^{M) be the finite codimension subspace 
of H^(M) defined by the condition 

f {{Vu,V(t)) + {u,(f>))dvM = V0Gker(£,/C). 

A Morrey-type argument by contradiction using (7.2) shows there is a constant 
C > such that for all u G H^{M), 

\ufdvM< \Cuf dvM + (p \J}Cuo\'^ dvY , (7.4) 
Jm Jm Jy 

where J = (1 + |^|)^/^. Now suppose {uk}^ C H^{M) is such that Cuk = fk ^ 
f G L2(M) and IC{uk)o = Sfc ^ cr G Pi7y^(F). (Note that by the definition 
(6.18) of /C, the range of is a subspace of PH^^iY)). Since the kernel is 
finite dimensional we may normalise Uk G H^{M), and then (7.4) and (7.2) 
show that {uk}i^ is bounded in H^{M). It then follows as above that there is a 
subsequence converging strongly in H^{M) to u, and that Cu = limfe— ►oo fk = f 
and ICuq = limfc^oo Sk = so the range of {C, /C) is closed. ■ 

The general boundary value problem 

ICuo = a on Y 

is solvable for u G H^{M) provided (/, cr) satisfies the condition (7.6) of the 
following main result. 
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Theorem 7.3 {f,a) e L'^{Y) x PH^^iY) lies in the range of {C,)C) (that is, 
(7.5) admits a solution u G H^{M)), if and only if 

[ if, 4>) dvM + / (cr, v(t>o) dvY = Q y (pe kev{£^,IC^) . (7.6) 

Proof: If n G H^{M) satisfies jCu = f and Kuq = a then u is also a weak 
solution. Condition (7.6) then follows directly from the definition 6.3 of weak 
solution, hence (7.6) is a necessary condition for solvability. 

To establish the converse, consider first the case a = 0. Thus we suppose 
/ G L'^{M) satisfies Jj^^if, 4>) dvM = for all 4> G ker(>C^, /C^), and we must find 
u G H^{M) satisfying Cu = /, K,uq = 0. 

By Lemma 5.1 the trace map ry '■ u i-^ uq is bounded, hence 

Hi := {u G H^{M) : Kuq = 0, and 

/ {{Vu, V(p) + (n, (j))) dvM = V G ker(£, /C)} (7.7) 

JM 

is a closed subspace of H^{M). The argument of Theorem 7.2 (ii) shows there 
is a constant C such that 

/ i\Vuf + \uf)dvM<cf \CufdvM (7.8) 

JM JM 

for all u G Hj(^{M). In particular, Jj^ \j0.u\'^ dvM is strictly coercive on Hf^, so 
the Lax-Milgram lemma gives u G i?^ satisfying 

/ (/, Ccj)) dvM = / (-Cm, £(p) dvM 

JM JM 

for all ^ G i?^. This equality also holds if G kev{C, /C), so * = Cu — f satisfies 
/ (*,>C(^) dvM = V G H^{M) , /C(^o = . (7.9) 

JM 

Lemma 5.8 and the identity 

/ £0) dvM = I {C)^, 0) dvM - i (j^*o, h) dvM (7.10) 

JM JM JY 

show that (7.9) is the weak form of the adjoint problem 

= 0, /C^*o = . (7.11) 
By (3.7), jC^ is elliptic with boundary representation 

= -Hdx + A + B), 

where A = -u^^Au since = A, and B = -v'^B^iy. By (6.8,6.9) the lead- 
ing terms in A are a*9j so A is elliptic on Y, and self-adjoint by (6.7). Since 
A{h'~^(f)ct) = —^a^~^(t>a if A(l)a = \a(f^a, we See that A satisfies the spectral 
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conditions, and spec A = — specA. (Note that in the usual case of Dirac opera- 
tors, A = A and the spectrum is symmetric). Now P := 1 — u^^Pv is a positive 
eigenspacc projector for A, with eigenvalues — for a G A~ U (A°\A), and the 
boundary operator satisfies 

K'^^q = {P + v-^K^v{1-P))-^o. (7.12) 

Since K = — z/~^i^V maps negative eigenvectors (of .4) to positive eigenvectors, 
it follows that 10^q = Q is an elliptic boundary condition for The bound- 
edness conditions (5.55,5.56) for K follow from the corresponding conditions 
for K. 

Since (£^,/C^) is elliptic and satisfies the conditions for Theorem 6.4, we 
conclude that ^ G H^{M) and ^ satisfies the strong form (7.11). 
Since ^ G ker(>C''', /C'''), assumption (7.6) with cr = gives 

/ {f,^)dVM = Q. 

Jm 

By construction ICuq = 0, so we may use « as a test function in the weak form 
(7.9) of the equation satisfied by giving 

/ {Cu, dvM = 0. 
Jm 

It follows from ^ = Lu — f that = and thus u is the required solution. 

Now consider the case a ^ 0. By Lemma 5.1 there is an extension v = 
ey((T) G H^{M) supported in a neighbourhood of y such that uo = cr, ||'y||//i(M) ^ 
2 II (7 II ^1/2. Let f = f — Cv and consider the equation 

Cu = /, /C«o = . (7.13) 
The previous case shows there is a solution provided / satisfies 

{f,tp)dvM = VVeker(£t,/ct) ci?^(M). 
Now (7.10) shows that for all G ker{C^,}C^), 

/ {f,ip)dvM = / {f,tlj)dvM- {v,d^tl:)dvM + f {vQ,yiiQ)dvY 

JM JM JM JY 

= / (/, V') dVM + f (o", i^V'o) dVY ■ 
Jm Jy 

Thus if (7.6) is satisfied then there exists a solution u of (7.13), and then 
u = u + V is the required full solution. This establishes sufficiency for the 
condition (7.6). 

■ 

We note two important consequences of Theorems 7.2, 7.3. 

Corollary 7.4 (7.5) admits a solution for all {f,a) G L^{M) x Pi?y^(y) if 
and only ifkev{C\}C^) = {0}. 



/ 

JM 
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Corollary 7.5 (>C,/C) : H\M) L'^{M) x PhI'\y) is Predholm. 



Proof: The argument of Theorem 7.3 shows that (£^,/C^) is eUiptic and thus 
has finite dimensional kernel by Theorem 7.2. Now (7.6) shows that the range 
of (>C,/C) has finite codimension. ■ 



8 Fredholm properties on complete noncompact man- 
ifolds 

In this section we establish conditions under which the Fredholm and ex- 
istence results of the previous section for the operator (>C,/C) : H^{M) 

may be extended to non-compact manifolds. This includes 
in particular, a generalisation of the solvability criterion (7.6) of Theorem 7.3. 
Results of this type may be applied to establish positive mass results in general 
relativity, for example. 

The non-compactness of M causes some difficulties not found in the compact 
case. A classical result [45, 65] shows that a Dirac operator D on a non-compact 
manifold is essentially self-adjoint on L?'{M). However, this elegant result is 
useless for our purposes, since it implies only that {{(p^Vcj)) : cj) G domP C 
L?{M)} is closed in the graph topology on L'^{M) x L?{M). This is weaker than 
the closed range property, which is necessary for useful solvability criteria. In 
fact, because L^(M) often does not encompass natural decay rates of solutions, 
the self-adjoint closure may not have closed range. In such cases the Dirac 
operator defined on (M) will not be semi-Predholm. This is shown explicitly 
in the following example. 

Consider the self-adjoint closure V : domV C L^(M^) L'^{M?) of the 
constant coefficient Dirac operator V = 7*0j and let / = Vu^ u = {\ — x)\x\~^i^, 
where x(r) is a smooth compactly supported function identically one around 
and ■0 is a constant spinor on R^. Clearly / G L^(IR^) but u ^ L^(R^), so in 
particular, u domD. However, / still lies in the closure of the range of P, 
since V{xru) = Xnf + Dxru ^ / in L^{R^), where xr{x) = x{x/R), but xru 
can not converge in L^(M^). Clearly {u,f ) graphP since u ^ L^(R^), and it 
can be shown (using the corresponding Schrodingcr-Lichncrowicz identity) that 
there is no u € L^(R^) satisfying Vu = /. Thus the self-adjoint closure V does 
not have closed range. 

In order to obtain an operator with closed range, it is thus necessary to 
enlarge the domain, which raises the question of determining the appropriate 
decay rate. We sidestep this problem by using the size of the covariant 
derivative as a norm. To obtain sufficient control on the L'^^^ behaviour, we 
then must postulate a weighted Poincare inequality (8.3). The existence of 
such inequalities can be established for the applications of most interest in 
general relativity; see Proposition 8.3 and §9. 

The elliptic estimate (7.2) plays a central role in the analysis over a compact 
manifold, but its noncompact analogue cannot be obtained directly by similar 
localisation arguments. However, in cases of geometric interest an identity of 
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Schrodinger-Lichnerowicz form (generalising (2.5)) is available, and can be used 
to construct suitable global estimates. 

The weighted Poincare and Schrodinger-Lichnerowicz estimates are the two 
additional ingredients needed for establishing solvability and Predholm proper- 
ties on a non-compact manifold. 

For ease of further reference, let us summarize the hypotheses which will be 
made throughout this section: 

Hypotheses 8.1 M is a non-compact manifold with compact boundary Y , 
which is complete with respect to a C°° background metric g. The case y = 
is admitted. The operator C = a^dj + b satisfies the global uniform ellipticity 
and boundedness condition 



for some rj > 0, for all V € Ex and all x € M . The coefficients of C satisfy the 
interior regularity conditions (3.4), and the boundary regularity and structure 
conditions of ^6, namely (6.2), (6.5), (6.7), (6.9). Let A be the boundary oper- 
ator and P its associated positive spectrum projection, as in ^6. The boundary 
operator K : {1 - P)L'^{Y) PL'^{Y) satisfies (5.55,5.56), and K. ,/Ct are 
defined in (6.18,6.20). The connection 



satisfies (3.18,3.19) and we note again that V need not be compatible with the 
metric on E — this is important in some applications. 

We may express C in terms of V by 



where (3 G L[^^(Af). Additional, rather weak, decay conditions will be imposed 
on (3 (8.8), on the negative part of the curvature endomorphism p (8.13), and 

on = 5(r + *r) in §9. 

Definition 8.2 The covariant derivative V on E over M admits a weighted 

Poincare inequality if there is a weight function w G L]Q^(Af) with ess infQW > 
for all relatively compact Q d M, such that for all u G C^(M) we have 



Here the length |Vnp is measured by the metric on E and the background 
Riemannian metric g on M, and dvM is the volume measure oi g. It is clear 
that the weight function w can be chosen to be smooth. 

The semi-norm 



(8.1) 



V = d-T 



(8.2) 



C = a^Vj + {h + a^Vj) = a^Vj + /3, 




(8.3) 




(8.4) 



on C^{M) may be completed to form the space 



H := II • lln-completion of C^r{E) , 



(8.5) 
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which consists of equivalence classes of M-convergent sequences in C^{M). 
The weighted Poincare inequality (8.3) ensures that an H-convergent sequence 
converges locally in L^, so the equivalence classes may be identified with cross- 
sections in the usual Lebesgue sense: with cross-sections having coefficient func- 
tions agreeing dfM"^.^ 

If there is no weighted Poincare inequality, then it may be that EI can not be 
identified with a space of Lebesgue-measurable cross-sections in this sense. For 
example, the trivial spinor bundle over M = x R with the flat connection 
Vj = di admits a global parallel spinor Viip = which is approximated in 
the M seminorm by -0^. = x(.x/fc)^j for x ^ C'^(M), x = ^ on [—1,1]. Now 
Im I^V'fcpc^^^M 0, but limfc^ooV'A; = V' / 0) so the H-equivalence class [0] 
contains V' 7^ everywhere. In other words, (8.4) does not define a norm on 
spinors in this example. This shows, inter alia, that (8.3) will not hold in all 
cases. 

More generally, a weighted Poincare inequality fails for manifolds of the form 
iV X R, where N is compact and itself admits a parallel spinor. It follows from 
the proof of Theorem 9.3 below that in such cases the orthogonal complement 
in EI of the subspace of all parallel spinors will admit a weighted Poincare 
inequality. Note also that the presence of a weighted Poincare inequality (8.3) 
does not imply there are no global parallel spinors — provides a simple 
counterexample . 

However, weighted Poincare inequalities can be demonstrated in many cases 
of interest. In the next section we will prove: 

Proposition 8.3 A covariant derivative V on E admits a weighted Poincare 
inequality if any one of the following conditions holds: 

1. there is a relatively compact domain CI G M and a constant c > such 



for allueC^{M); 

2. there are no nontrivial globally parallel sections (Vu = ^ u = 0); 

3. M has a weakly asymptotically flat end M (see Definition 9.4), with 
dimM > 3; 

4- M has a weakly asymptotically hyperholoidal end (see Definition 9.9), with 
dimM > 2. 

When M is non-compact, the global Garding inequality (generalizing (7.2)) 
cannot be constructed from local estimates. Motivated by some classical and 
fundamental identities, we instead introduce the following definition. 

Definition 8.4 The operator pair {£,,}€) admits a Schrodinger-Lichnerowicz 

estima( (- if III ere is C > and a non-negative function p such tha^ for all 

**Thc function p here should not be confused with the energy density arising in general 
relativity: in Section lip will be zero. 



that 




(8.6) 
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C-^ \VufdvM< / {\jC.u\'^ + pIu]"^) dvM + <p \J}Cuo\^ dvY , (8.7) 
Jm Jm Jy 

where J = (1 + \A\)'^/'^ . 

Lemma 8.5 Suppose that the Schrddinger-Lichnerowicz estimate (8.7) holds for 
all u G C^{M) with p and P = C — a^Vj satisfying p G L^oc"^ ' ^ -^loc '^'^^ 

lim sup P±\^ < oo , (8.8) 

where {Mjiju^oo is an exhaustion of M. Then C : M ^ L'^{M) is bounded and 
(8.7) holds for allu&U. 

Proof: It will suffice to show that the individual terms of the right-hand-side 
of (8.7) are bounded by Now 

f \£uf<c[ |Vu|2 + 2/ (8.9) 
Jm Jm Jm 

and we use (8.8) and (8.3) to estimate 

/ {p+m')\u\' < sup p±}^ f 

JM\Ma M\Ma W J M\Ma 

P + \(i? f 

< sup f 

M\Mr W JM\Mji 

< C [ \Vuf , 

Jm 



for some R < oo. Let XR ^ C^{M) be a cut-off function with support contained 
in M2R, Xi? = 1 on Mr. Then 

ip+\f3f)\nf < [ ip + \Pf)\xRuf 

Mr JM2R 

< (l|p|lL«V2(M2fl) + \\I^\\l"*{M2r)) WXRuWln* ^M2r) ■ 



Applying the Sobolev inequality for V on the compact set M2R and the weighted 
Poincare inequality show that the last term is controlled by Jj^ |Vnp. Finally, 
the iC-bound (5.55) and the restriction Lemma 5.1 show that the boundary 
term is also controlled by /^^ |Vnp. ■ 

Schrodinger-Lichnerowicz identities hold for many common examples, and 
can easily be adapted to produce estimates of the form (8.7). We will not 
attempt to give general conditions which imply such inequalities — it is simpler 
to ask only that (8.7) be established separately in any particular case of interest. 
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For example, consider the classical Dirac operator P of the metric g as in 
§2, on a non-compact spin manifold M. Combining (2.6) and (2.8) gives 

[ IWI' dvM = [ i\V^f - jRigM^) dvM + i {^Po, {Vy + '^Hy)^Po) dvy , 
Jm Jm Jy 

(8.10) 

for any spinor field on M. Suppose the boundary operator is /C = P+, 
the orthogonal projection onto the positive spcctrTim eigenspinors of Dy. If 
the boundary mean curvature Hy satisfies Hy < iyi67r/Area(y), then the 
argument in §2 shows that the boundary term in (8.10) is not greater than 

^^{P+iPo,VyP+'iI)o) dvy < II^V'o||J^i/2(y) ' 

and (8.7) follows immediately, with 

p = max(0,-ii2(5)) . (8.11) 

Since /? = in this example, the inequality holds for all u G EI provided p 
satisfies (8.8). For general mean curvatures Hy G L°°{Y), note again that 



<* {'ll^o,Vytjjo) dvy < ||P+Vo||li 
Jy ^* 



/2(y) • 

If Hy G L°°(y) then §y Hy\i^o\^ < ||//y |Uoo(y) || J^oll^i/a^^^- Using a frac- 
tional Sobolev inequality, the control on Hy may be weakened to Hy G L'f{Y), 
p = n—1 for n> 3 and p > 1 for n = 2. Lemma 5.1 shows that ||«/V'o||j|^i/2^y-| < 

c||V'o||_ffi(yx/')) where ip = xV' ^^'^ X = x{^) is a cutoff function supported in 
/' = [0,(5/2], as in the proof of Theorem 6.4. Now Corollary 6.5 shows that 

C'HfHUYxI') < f d^V'l' + IV'I') dVM + i I JP+V'ol' dvy, 

which provides the required Schrodinger-Lichnerowicz estimate (8.7). 

In applications, a Schrodinger-Lichnerowicz estimate is usually obtained in 
the special case of homogeneous boundary data (/Cmq = 0). The above trick 
shows that the estimate in the homogeneous case implies the general estimate 
(8.7): 

Lemma 8.6 Under the hypotheses of Lemma 8.5, suppose there is C > such 
that for all u with ICuq = we have 

C-^ [ \VufdvM< f {\jCuf + p\uf) dvM , (8.12) 
Jm Jm 

for some p. Then there is C > such that (8.7) holds for all u EM. 

Proof: Suppose n G H and let u = u — xu, where x = x(^) ^ C°°{M) is 
a cutoff function supported in F x /' as in the proof of Theorem 6.4. Then 
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JCuq = so (8.12) applies to u, giving 

/ \Vu\''dvM < 2/ {\Vu\^ + \V{xu)\^)dvM 

JM JM 

< c[ {\Cu\'^ + p\uf) dvM + 2 f \V{xu)fdvM 

JM JyxI' 

< c[{\Cuf + {p+\dxf)\uf)dvM + 2[ \V{xu)fdvM . 

JM JYxI' 

Now it follows easily from Corollary 6.5 that 

f W{xu)?dvM< I {\Cu\^ + \u\^)dvM+ i \J^UQ\^dvY, 

JYxI' JyxI JY 

which gives the required inequality. ■ 

Theorem 8.7 Under the hypotheses 8.1, suppose (M, V, C, K.) admits a weighted 

Poincare inequality (8.3) and a Schrddinger-Lichnerowicz inequality (8.7) with 
p and (5 satisfying the conditions of Lemma 8.5. If p & Lf^^{M) for some 
p > n*/2, and if 

lim sup - = , (8.13) 

R^oo M\Ma W 

where {Mii^r-^^, is any exhaustion of M, then 

{£, /C) : M ^ l2(M) X H^^iY) (8.14) 

is semi-Fredholm. 

Proof: Lemma 8.5 gives £u G L'^{M) for u G H. We first show the unit ball 
in the kernel is compact. Let {uk}'^^^ be a sequence in the kernel of {C,IC), 
normalised by ||iifc||H = 1- Weak compactness of bounded sets in M shows there 
is n G EI and a subsequence, which we also denote by Uk, such that Ufc — ^ n G H 
and \\u\\m < liminf ||nfc||H = 1- 

Since (8.13) is independent of the choice of exhaustion, we may suppose 
for definiteness that Mr = {x e M : d{x) < R} where d{x) is the smoothed 
distance function from some fixed base point. Let x € satisfy x(x) = 1 

for X < I, x{^) = for x > 2 and < x{^) ^ 1) lx'(2^)l ^ 2 for all x. Then 
the functions xr{^) = x{d{x)/R) form support functions for the exhaustion 
Mr which satisfy suppxi? C M2R, XR = 1 on Mr and \dxR\ < 2. Using the 
weighted Poincare inequality we have 

/ \'^{XRUk)\'^ dvM < 2/ \dxR\'^\uk\'^ dvM + 2 {Vukl"^ dvM 

JM JM2r\Mii J M2R 

< 2(1 + 2 sup w~'^) / {Vukl'^dvM, 

M2r\Mr JM 

which shows that for any R > 1 the sequence XR'^k is bounded in H^{M2r). 
Since XR^k XR.^ in H^{M2r), the Rellich lemma implies XR'i^-k ~^ XR^ 
strongly in L1{M2r) for any g < 2 = 2n/(n — 2) and any R> 1. 
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Applying (8.7) to any difference Uj — Uk gives 

/ \V{Uj - Uk)\'^ dVM < / p\Uj - Ukl"^ dVM 
JM JM 

< \\p\\LP{Mii)\\Uj - '"fellL(MH) 

+ sup — / \uj — ui;\'^ w dvM (8.15) 

M\Mh JM 

wfiere, since p > n*/2, we liave q = 2p/{p — 1) < 2. Now (8.3) and llufeHn = 1 
combine to sliow tliat 



/ I "J 

JM 



Uj — Uk\'^ wdvM < 4, 



so by (8.13), for any e > there is i? = R{e) such that the second term of (8.15) 

is less than e/2 for all j, k. Since Uk converges in L'^(Mji) there is N = N(e, R) 
such that the first term is less than e/2 for all j,k > N. This shows Uk is a 
Cauchy sequence, hence strongly convergent to u, in H. 
As noted above, ||>C'u||^2(jv^) < CHuHh and thus 



/ \£u\'^dvM = / \jC{u - Uk)]"^ dvM 

JM JM 



< C \V{u ~ Uk)\'^ dvM 

JM 

— > as — >■ OO , 

which shows that Cu = 0. Similarly, since /C : ijy^(F) ^ ijy^(y) is bounded, 
for any uEMwe have 

^ \ JJCuo\^dvY < c\\JCuo\\j^i/2^y^ < cA;||no||^i/2^^^ 
< C [ \VufdvM , 

JM 

by (5.55) and the trace lemma 5.1. Choosing u = u — Uk gives 
® \JKuq^ dvY < c I \V{u — Uk)^ dvM = o{l), 

JY JM 

which shows also that KLuq = 0. Thus u G ker(>C,/C) and the kernel is finite 
dimensional. 

To show the closed range property, observe that by (8.13) and (8.3), the 
elliptic estimate (8.7) may be strengthened to 

I (|Vt(|^ + |u|^u)) dvM < / \Cu\'^dvM+ / p\u\'^dvM+f \JICuo\'^ dvy , 

JM JM Jn JY 

(8.16) 

for some relatively compact domain CI (e M. Now we claim there is a constant 
C > such that 

/ p\u\'^dvM <c( I \jCu\'^dvM + i \ JlCuQ\^dvY] , (8.17) 

JCl \JM JY J 
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for all u G H such that 

/ (Vu, V(/)) dvM = V G ker(£, /C) . (8.18) 
Jm 

Suppose (8.17) fails, so there is a sequence Ufe G H, = 1, 2, . . ., such that (8.18) 
holds for each Uk, and 

/ /9|ufepdt;M = 1, / \jCuk\'^dvM + f \JIC{uk)o\'^ dvy < l/k. 

Jci JM JY 

The sequence is bounded in H by (8.16), so by passing to a subsequence we may 
assume Uk converges weakly to u G H and strongly in L^(n), q = 2p/{p—l) < 2 
as before. Applying (8.16) to Uj — shows the sequence is Cauchy and thus 
converges strongly in H. It follows that 

/ \£u\'^dvM + f \JICuo\^ dvy = 0, 

JM JY 

so n G ker(/^,/C). Strong convergence shows that (8.18) is also satisfied by u, so 
testing (8.18) for u with (j) = u shows that u = However, strong convergence 
in L'^{Q) shows that J^p\u\^dvM = 1, which is a contradiction and establishes 
the claim (8.17). 

Combining (8.17) with (8.16) gives 

/ (|Vtxp + \u\'^w) dvM <C [ ICul'^dvM + f | J/C«oP dvy] (8.19) 

JM \JM JY J 

for all -u G H satisfying (8.18). Now suppose tXfe G H is a sequence such that 
Cuk = fk ^ f ^ LP'{M) and IC{uk)o = Sfc — >■ s G hI^^{Y). These convergence 
properties are retained if wc replace by + yk for any convergent sequence 
Uk G ker(Zl,/C), so we may assume the Uk all satisfy (8.18). In particular, 
applying (8.19) to uj — Uk shows that Uk is Cauchy in H and converges to u 
satisfying jCu = /, JCuq = s. This shows (>C,/C) has closed range. ■ 

By Definition 6.3, u is a weak solution of 

Cu = /, /Cuo = C7 , (8.20) 

for f eL'^{M),ae PH^^iY), if n G Ll^{M) and 

/ {u, dvM = / (/, ^) dvM + <p {(T, y4>o) dvy , (8.21) 

JM JM JY 

for all 4> G Hl{M) such that = 0. Similarly, the argument of Theorem 7.3 
shows that the weak form of the adjoint problem 

jC^'u = g, /C^uo = T , (8.22) 

for g G L2(M), r G PHI^^{Y), P=1- v-'^Pu, is that u G L\^^{M) and 

/ {u,£(j))dvM= {g,(t)) dvM - f {T,iy~^(j)o) dvy (8.23) 

JM JM JY 

for all ^ G i^c(M) such that /C0o = 0. 

We now extend the solvability criterion (Predholm alternative) of Theorem 
7.3 to the non-compact case. 
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Theorem 8.8 Under the conditions of Theorem 8.7, suppose the formal adjoint 
{C'^,)C'^) also satisfies a Schrddinger-Lichnerowicz estimate (8.7) with the same 
covariant derivative V and with a curvature term p satisfying (8.13). Then the 
system (8.20) with {f,a) G L^{M) x PH^/'^(Y) has a solution u e M if and 
only if (/, a) satisfies 

[ if, 4>) dvM + i {(T, lycPo) dvY = 0, (8.24) 

for all 4> ^Mn L^{M) satisfying C'^cp = 0? = 0. In particular, the system 

(8.20) is solvable for all (/, a) G L'^{M) x PH^/^{Y) if and only if there are no 
^ eMn L'^{M) satisfying C^-^ = 0, K.'f-^o = 0. 

Remark 8.9 We emphasise that in Theorem 8.8 it is not necessary to impose 
conditions on p other than (8.13), and no conditions on the map $ := — *a*Vi 
are needed. 

Remark 8.10 See Theorem 11.9 for an example where £^ / £, with jC^ sat- 
isfying two Lichnerowicz-Schrodinger identities with respect to two different 
connections. 

Proof: The necessity of (8.24) follows immediately from the weak form (8.21). 
To show sufficiency, the argument of Theorem 7.3 applies to reduce to the case 
(7 = 0, which we now consider. 

Let Mjc = {u£M: Kuq = 0}. The elliptic estimate (8.7) gives 

/ \Vu\^dvM<c! {\Cu\^ + p\u\^)dvM , VuGMx;. 

JM JM 

The arguments used to show (8.16) and (8.17) apply and give 

/ |V«p dvM <C {Cul"^ dvM V It G Hyc, (8.25) 

JM JM 

where we define 

:= {« G M/c : / (Vu, Vcfy dvM = V ^ G ker(/:, JC)} . (8.26) 

JM 

Thus the bilinear form u ^ \Cu\^ dvM is strictly coercive on the Hilbert 
space H^Cj and for each / G L?{M) the map (j) — > Jj^if, C(t>) dvM is bounded on 
Ha:- The Lax-Milgram lemma shows there is n G Ha: satisfying 

/ {Cu, £(j)) dvM = / (/, dvM y (peMjc ■ 

JM JM 

Thus setting = Lu — f we have 

/ {^,C(P)dvM = V0gHk, (8.27) 

JM 
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since (p £ ker(£,/C) will also satisfy the relation (8.27). Lemma 8.5 shows that 
* G L2(M) and from Definition 6.3 and (8.27) we see that is a weak solution 
of 

= 0, = . 

If there are no such non-trivial then Cu = /, and u is the required solution. 
The arguments of Theorem 7.3 show that (jC^K.^) is elliptic and Theorem 6.4 

applies to show ^ G H^^^(M). Let Mji be the exhaustion of M constructed 
in Theorem 8.7, with associated cutoff functions XR ^ C^{M), and let ^'^ = 
Xk^ G HI{M) C H. The assumed Schrodinger-Lichnerowicz estimate (8.7) for 
{d^ ,K.^) gives (with curvature term p) 

f |V(*fe - dvM<C [ (\d{^k - + P\^k - dvM . (8.28) 

Since = we have 

/ \L\^k-'^ltdvM<c\ (|cixfe|'l*l' + |c?Xi|'l*l') c^^^M^O, 

JM JM 

because * e L'^iM), \dxk\ < 2 and suppdxfc C M2k\Mk. Now 
/ p\^k-^l? dvM <e i \^k-^l\^wdvM 

JM JM 

by the condition (8.13) on p, for sufficiently large fc, I. By the weighted Poincare 
inequality (8.3), this is in turn bounded by e times the left side of (8.28) and 
may therefore be discarded in (8.28) by choosing e sufficiently small. It follows 
that ^'fc is a Cauchy sequence in H, so * G ElnL^ and thus d^ = 0, K)^q = 0. 
If there is no such 7^ then Cu = /, and u is the required solution. More 
generally we have Cu = f + , u £ H/c, and since Jj^{Cu, ^) dvM = by (8.27), 
the condition (8.24) (with a = and (p = ^) shows that ^' = and we have 
solved Cu = /, as required. ■ 



9 Weighted Poincare Inequalities 

Before proceeding with the analysis, define the symmetric part of the con- 
nection V by the formula 

(0, r^{X)^) := ^ [X{cp, ^) - Vx^) - (Vx0, V')) , (9.1) 

for all smooth sections (pjip oi E and all smooth vector fields X: One easily 
checks that (9.1) defines a linear map r^{X) from fibers of E to fibers of E, 
symmetric with respect to the scalar product (•,•), with the map X r^{X) 
being linear as well. Clearly, V is compatible with (•, •) if and only if L*^ vanishes. 
If r is defined by (8.2), then 

= i(r + *r) . 

We establish Proposition 8.3 via a special case, based on an argument of 
Geroch-Perng [32]: 
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Lemma 9.1 Let be any two relatively compact domains in M, and assume 
that 

G UC{M) . (9.2) 
There is a constant e > such that for all sections u € Hl^^{M) of E we have 

e l\u\^dvM^ I \u\'^ dvM + / \Vu\'^dvM- (9-3) 
Jo. Jn Jm 

Proof: Let q be any point of f2, fix p G O and let be small enough that the 

5f-geodesic ball B{p, rp) of radius and centred at p, lies within VI. Let X be a 
compactly supported vector field, such that the associated flow (pt satisfies 
(j)i{B{p,rp)) D B{q,rq) for some > 0. (Since M is C°° and connected, it is 
always possible to construct such an X.) Let fit = (f)t{B{p,rp)). 

By direct calculation and Holder's inequality we have, for any u G Hl^^{M), 



d_ 
di 



/ \u\'^dvM = / {2{u,{Vx +T^)u) + \u\'^diVgX) dvM 



< c 



where C depends on ||divgX||ioo . By the Sobolcv inequality in the 

coordinate ball Vlt for functions, ||/||j^2(s^^) < C{\\df\\L2(^^) + ||/||L2(!^i))- Ap- 
plying this to / = \u\ gives ||^^||j^2(f^^) < C'(||-C>n||i2(n,) + \\u\\l2(^^)), where D is 
any metric-compatible connection. Since G L"* may be written as Fi + 
Fi G L°°, ||F2||^n* < e, the Sobolev inequality gives 

11^^11^2(0,) < C (||V?x||i2(no + hllL2(nt)) ' 
for some constant C depending on F. Defining F(t) = J^^ \u\^dv^, we have 

;F{t) < CF{t) + C [ \Vu\'^dvM , 



d_ 
dt' 



and Gronwall's lemma gives F{1) < e'-^{F{0) + \Vu\'^dvM)- Thus there is 
e > such that 

e / {ul"^ dvM < / {ul"^ dvM + / \Vu\'^ dvM ■ 

JB{q,rq) Jq Jm 

Since Cl has compact closure, it is covered by finitely many such balls B{q,rq) 
and (9.3) follows. ■ 

Corollary 9.2 Under condition (9.2), if there is a domain c M and a 
constant e > such that 

el \u\^dvM< / \Vu\^dvM (9.4) 
Jn Jm 

for all u G C^(M), then M admits a weighted Poincare inequality (8.3). 
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Proof: By paracompactness and Lemma 9.1, there is a countable locally finite 
covering of M by domains $1^ and constants 1 > > 0, A; G Z+, such that for 
each k, 

€k / {ul"^ dvM < / {ul"^ dvM + / {Vul"^ dvM ■ 
This is in turn bounded uniformly by (9.4), so the function 

^(^)= E ^ (9-5) 

is bounded, strictly positive, and satisfies 

/ \u\'^wdvM< / |VupdvM, 
Jm Jm 

which is the required weighted Poincare inequality. ■ 

This establishes part (i) of Proposition 8.3, and we next turn to the proof 
of part (ii). 

Theorem 9.3 Suppose that M has a locally finite cover such that 

Vi = di-ri, with FiGLjJ,;. (9.6) 

If there are no global V -parallel sections of the bundle E, then M admits a 
weighted Poincare inequality. Equivalently, if M does not admit a weighted 
Poincare inequality then M admits a global V -parallel section. 

Proof: Assume M does not admit a weighted Poincare inequality, so by Corol- 
lary 9.2, for each domain $7 (£ M and each constant e > 0, there is n G i?j^^(M) 
such that (9.4) fails. In particular, fixing for each A; > there is € H^^^{M) 
such that 

l\uk\^dvM = l, I \Vuk\^ dvM < k-^ ■ (9.7) 
Jn Jm 

It follows that V-Ufe — strongly in L'^{M). Under (9.6) Rellich's lemma holds, 
so there is a subsequence converging strongly to n G L'^(ft). Then Vu = and 
u 7^ in O. 

Now let Mj, J = 1, 2, ... be the exhaustion of M from Theorem 8.7, and let 
Uj G H^{Mj) be the corresponding parallel spinors, constructed in the preceding 
paragraph. Since Uj ^ there is Mj d M^- such that ^^j^, \uj\^ ^ 0. Lemma 9.1 

applied with Mj replacing M shows there is r]j > such that for all v G 

rjj I \v'^ dvM ^ I \v'^ dvM + / \Vv'\^ dvM ■ 
Jm', Jmo Jma 

In particular this implies | Uj p dvu 7^ and we may impose the normalisa- 
tion dvM = 1- By Rellich's lemma there is u\ G H^{Mi) and a subse- 
quence, also denoted by uj, such that Uj — ui in H^{Mi) and = 1, 
Vui = 0. 
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Again by Lemma 9.1, for each k>l there is > such that 

Cfe / l^l^ dvM < / bp dvM + / iV'up dvM , y V e Hl^^{Mk+i) . 

JMk JMi JMk+i 

Setting V = Ui — Uj, i, j > k, shows that the sequence Uj is Cauchy in L?'{Mk) 
and therefore converges strongly in L'^{Mk) for all A; > 1 to some nontrivial 
u e i^iocW, and Vtx = 0. ■ 

Another application of Corollary 9.2 leads to Proposition 8.3 part 3, for 
asymptotically fiat manifolds. In fact the proof works for a much broader class 
of manifolds: 

Definition 9.4 A weakly asymptotically flat end M <Z M of a Riemannian 
manifold M with metric g is a connected component of M\K for some compact 
set K, such that M ~ M"\i3(0, 1) and there is a constant rj > such that 

for all X e W\B{0, 1) and all vectors ^ G W. 

Theorem 9.5 Suppose {M,g) is a (connected) Riemannian manifold of dimen- 
sion n > 3, g E C^{M), and M has a weakly asymptotically flat end M. Sup- 
pose also the connection Vi = di — Fj on E satisfies F G Lf^^{M) and the decay 
conditions 

where is the symm,etric, scalar product incompatible, component ofV defined 
by Equation (9.1). Then M admits a weighted Poincare inequality. 

Remark 9.6 The restriction dimM > 3 is rather harmless as far as the ap- 
plications to the positive mass theorems are concerned, since the notion of 
asymptotic flatness for two dimensional manifolds, relevant to general relativis- 
tic applications, has to be defined in a completely different way. An adequate 
analogue of mass here when dimM = 2 is provided by the Shiohama theo- 
rem [57]. 

Remark 9.7 The decay condition (9.8) is independent of the choice of flat 

background metric gij = Sij: Equation (9.1) shows that is a tensor. By 
comparison with the ^'-distance function from any chosen point p, the function 
r is equivalent to this distance function, which implies the result. 

Remark 9.8 The proof below establishes the inequality (9.4) for spinors sup- 
ported in O := \ -6(0, R) for some R without assuming that F G U^*^{M). 

Proof: Let r = (E(a;')^)^^^ e C°°(M) and x = x{r) e C^{M) satisfy, for 
some i?o > 1 and k > 10, 

Xir) = ^^, 2R,<r<ik-l)R, 
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and x(r) = 1 for r > kRQ, x(r) = for r < Rq. Then |x'(r)| < 2/(rlogA;), so 
for any section u G C^{M) 

[ \V{xu)\'^dvM f \Vu\'^dvM + jr^TTo f \\ufdvM. (9.9) 

Now Ao(r2-") = for r > 1 in M**, n > 3, so for any v G Cl{W^\B{Q,Ro)) we 



have 



= - / ai(ai(r2-")|t;|V"-2)dx 

= {n-2f r-'^\v\'^ dx + (n - 2) r^^2{v,iVr +T^)v) dx , 

where = r~^x*rf and lengths are measured by g and the metric on E. Using 
Holder's inequality we obtain 

~ / r-2|^;pda:< / | Vt;^ + (n - 2) / r''^\vf\Tf\dx. 

4 Jr" Jr" 7r" 

The Sobolev inequality in M", n > 3, 

l-2/n 



< 

< 4 



/ \vfdx) <Cs \Dvfdx , 

Jr" J JR" 

where D = 'V + is the metric-compatible connection, gives the estimate 

/ \Dvfdx < 2 /" {\Vvf + \vf\T^f)dx 
7r" Jr" 

2/ |Vt;|2da; + 2C5||r^||i„(R„\5(o,i?,)) / 
7r" 7r" 

JR" 

provided 2C5||r'^||in(jjnyB(o,ilo)) — 5- ^^ow (9.8) implies there is i?o < 00 such 
that this condition will be satisfied, so for any v G C^{W-\B{0, Ro)) we have 

/ r-^\rf\\v\'^dx < ||r-^r^||i„/2(u„)C5 / \Dv\'^ dx 

JR" JR" 

< 4C5||r-ir^||^„/2(Kn) / \Vv\^dx. 
Hence there is e > such that for all v G (M n {r > Rq}), 

e lr~'^\vf dvM < l\Vv\'^ dvM ■ (9.10) 
Combining (9.10) with v = xu and (9.9) gives 

/ r~^\u\^ dvM < lr~^\xufdvM 

J{r>kRo} Jm 

C Lwixu^dvM 

JM 



< C I \Vu\'^dVM+ ^,.0 I 

Jm (log ky Jm 



r dvM , 



'm 
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where now \Vu\^ = g^^ {V iU,V ju) . If k is chosen so that C/(logA;)2 < ^ then 
the last term may be absorbed into the left hand side, giving 

/ r~'^\u\^ dvM < C \Vu\^dvM- (9.11) 

Jr>kRo Jm 
Lemma 9.1 now applies and gives the required weighted Poincare inequality. ■ 
In order to prove part 4. of Proposition 8.3 the following Definition is needed: 

Definition 9.9 A weakly hyperboloidal end M G M is a connected component 
of M\K for some compact set K, such that M ~ (0, xq) x M , where {J\f,h) is 
a (boundaryless) compact Riemannian manifold with continuous metric h, with 
g\^ being uniformly equivalent to 

g = x~'^{dx^ + h) . 

Here x is the coordinate running along the (0, xq) factor of {0,xo) xAf. 

We have the following hyperboloidal counterpart of Theorem 9.5: 

Theorem 9.10 Suppose {M,g) is a (connected) Riemannian manifold of di- 
mension n > 2, (7 € C^{M), and M has a weakly hyperboloidal end M. Sup- 
pose also the connection Vi = di — Fj on E satisfies T G LJJj^(M) and the decay 
condition 

limsupl^rfl < ^^^^ (9.12) 

in M, where is the symmetric part ofVg^, with norm understood as that of 
an endomorphism of fibres of E. Then M admits a weighted Poincare inequal- 
ity. 

Proof: This is essentially McKean's inequality [49]; we follow the proof in [22]. 
Let, first, / be a function in C^([0,xo] x J\f) with / = at {x = 0}; we have 

f{x,v) = 2jj{s,vfJ^ds 

Here we use the symbol v to label points in A^. Integrating on [0, xq] x A^, 
a change of the order of integration in x and s together with some obvious 
manipulations gives 

/ f'^x'"^dxd^h < - — ^— ^ / (x^i x~"^dxdnh 

J[OM^-^ ~ V i[0,a;o]xAr \ OX J 

< , ^,,2 / g{df.df)x-''dxd^xh.{^.l^) 

This is the desired inequality on M with metric g for functions, with weight 
function w = {n — 1)^/4. The result for general weakly asymptotically hyper- 
boloidal metrics and for functions follows immediately from the above, using 
uniform equivalence of g with g on the asymptotic region, and using Lemma 9.1. 
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Let, finally, t; be a smooth compactly supported section of a Riemannian 
bundle with not-necessarily-compatible connection V. Let (p be any smooth 
compactly supported function equal to 1 on the support of v, set 

fe = (l)\/e+{v,v) . 

We have 

= K^|=(. + <„,„))+/<^i(Yf + r»(..(v. + r» 



dx 



e+ {v,v) 
< e|#p + 02|(V. + rf)t;p. 



The first line of (9.14) yields 
/ x~"'dx d/^h = / (/^{e + {v,v)) x~'^dx 

JM Jm 

^ r^U2 [ x^e\dcf>\^ + cp\V, + Tl)v\^) x-''dxdf,h. 
(n - ly Jm 

Passing with e to zero gives 

/ {v,v) x~'^dxdfj,h < 7 TT^ / x^|(Va; + rf )up x~"'da;d/x^ 

Jm — Ij Jm 

< , ^,xo / {{\^\)\^vf^ + x^{\^b)\Vlvfdx\x-'^dxdiih 
[n-iy Jm \ « / 

for any S > 0, and if condition (9.12) holds the last term can be carried over to 
the left hand side, leading to 

C"M {v,v)dvM < / \Vv\ldvM ■ 

JM JM 

Lemma 9.1 gives then the desired inequality, with a weight function w equal to 
1/C in the asymptotic region. ■ 



10 Examples and Applications 

The structure and regularity conditions may be readily verified in situations of 
interest, which we illustrate by considering the Dirac operator examples of §2. 

Suppose M is a Riemannian spin manifold. Fix local coordinates (x^) and 
a local orthonormal framing = e'^d^ of the tangent bundle TM, and let 
(f)i, I = 1, . . . ,dimS' denote an associated spinor frame, determined by some 
choice of representation c : C£n — End(S'). The Dirac operator V defined by 
(2.2)-(2.3) is 

^ = efya^ - i^.,(9^)e^S^7^' , (10.1) 

where the skew-symmetric matrices 7* = c{ei) G End(5) are constant in the 
local spinor frame and satisfy the Clifford relation 

+ y y = _25»i . (10.2) 
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Clearly V has the form a'^d^ + h where 



n (10.3) 

i,i,fc,/i=l 

If gnu S Wj^'" n (7[q^ then by the Gram-Schmidt construction, the local or- 
thonormal frame may be chosen so that the coefficients also satisfy ef € " ^ 
C|^p, and then a^,6 also satisfy (3.4), cf. Proposition A. 9, Appendix A. The 
Clifford identity shows l^^a'^Vl'^ = g'^''Cn^v\V\'^ pointwise, which implies (3.5) 
so Theorem 3.7 may be applied to establish the H^^^ interior regularity of 
weak solutions. 

For boundary regularity we assume there is a diffeomorphism of y x / with 
a neighbourhood of y x {0} = dM in M, such that for any chart of Y, the 
associated chart of adapted coordinates {y^,x), x >0 satisfies 



9xAiy,0) = 0, ^ = l,...,n-l 

fc(y,o) = 1, 



(10.4) 



and 

9^^u\y e (Y) n C°(y) . (10.5) 

(Recall that the symbol "|y" stands for "|yx{o}"-) 

For smooth metrics this follows easily using Gaussian coordinates about Y. 
Geodesic uniqueness may not be available in the more general case of a VF'^+^'P 
manifold with metric g^u G VF'^'^ for k > n/p; in this case the existence of 
boundary coordinates satisfying (10.4) is guaranteed by Proposition A.IO. 

Then (6.5) follows from the regularity conditions on g^^, (6.2) follows from 
(10.5), and (6.9) follows from (10.4), the Clifford relations and the skew-symmetry 
of the 7*. In terms of an adapted frame e^, where e„|y = —dxn is the outer 
normal at Y = dM, the boundary operator (6.10) may be taken as 

n— 1 n— 1 

A=vy = -Y, efTV^M + i E ^^.(9M)e^Vyy , (lo.e) 

i=l i,j,k=l 

by an appropriate choice of 6o in (6.10). Again note that other choices of A are 
possible, such as Vy + F for any function or symmetric endomorphism F. 
The boundary condition (2.13) 

P+i/j = a, ae P+hI''^{Y) , (10.7) 

where Pj^ is the orthogonal projection to the positive cigenspace of A, corre- 
sponds via (6.14) to i^T = 0, which clearly satisfies the conditions of Theorem 
6.4 for K,K^. Prom Theorem 6.4 it follows that the weak-strong property and 
the elliptic estimate (6.21) hold at the boundary for the boundary condition 
(10.7). 
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Boundary operators of the type (2.18) were used in [30, 31, 47], for example. 
More generally suppose there is an endomorphism e : S ^ S acting on sections 
of which satisfies 



e^ = l, *e = e, eA + Ae = 0. (10.8) 

Note that if there is a splitting E = Eq® Fq such that A has the form (4.10) 

then e = ^ will satisfy (10.8). Conversely, if (10.8) holds then setting 

Eq,Fq equal to the ±1 eigenspaces of e shows that A may be written in the 
supersymmetric form (4.10). 

If Aij) = \tp then ^(e'0) = —eAij: = —Xeip, so the spectrum of A is symmetric 
and e interchanges the positive and negative eigenspaces. Let P+ (P-) be 
orthogonal projection to the positive (negative) eigenspaces, and consider the 
eigenspace splitting L^{Y) = iJ_ Hq Hq, where Hq = kei A and 
Hq = Hq Hq is the decomposition into ±1 eigenspaces of e. There is an 
isometric isomorphism e : H^ — > iJ_ such that e has the block decomposition 



.-1 



, acting on = i7+ iJ_ 



It follows that the action of K.+ = i(l + e) is given by 



H+ 



H. 



• 



/C+V' = 5(l + e) 



p 




q 








r~ 





\{ep + q) 




(10.9) 



Hence the boundary condition /C+V' = cr where a = ea = e(y+, ctq", 0], is 

equivalent to the component conditions \{p + e~^q) = C7+ and 
may be expressed in the form (6.14) if we define P, K by 



This 



K 



(10.10) 



since then 



Pi) - K{i - p)ip 



p 




e ^q 




" 2o-+ " 

















+ 































(10.11) 



is equivalent to /C+V = cr as above. Clearly K is bounded on both L'^{Y) 

1/2 

and Hj {Y) as required by the regularity theorems of §5. Similarly we find 
that the boundary condition JC-ij) = cr is equivalent to \(!P — e~^q) = cr_ and 
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ctq , so an equivalent spectral projection condition may be constructed 



using P = P+ + and K = -K, 



Pip - K{1 - P)'4) 



so 



p — e 




' 2(7+ " 


















r~ 




. ^0 . 



(10.12) 



11 Positive mass theorems 

Our motivation for the present work arose from positive energy theorems, and 
we shall present here some such theorems which follow from our work elsewhere 
in the paper. The main point is to give a complete proof of Hcrzlich's inequality, 
cf. Theorem 11.7 below. In addition, our remaining results in this section im- 
prove the previous ones [8, 11, 19, 20, 25, 37, 38, 51, 54, 62, 64, 67] in two respects: 
1) the low differentiability of the metric; 2) we do not require M to have a com- 
pact interior. This second aspect of our results is critical for some applications 
of the positive mass theorem to black holes uniqueness theory [17,23,24,53]. 
Some of the arguments already presented in other sections will be repeated in 
the proofs below, whenever useful for the clarity of presentation. 

A set (M, g, K) will be called initial data for Einstein equation if (M, g) is a 
three^ dimensional Riemannian manifold, and K is a, symmetric tensor on M; 
this is a slight abuse of terminology as we are not requiring any constraints to 
be satisfied. Only complete (M, 5)'s will be considered, with boundary either 
compact or empty. Given such a triple we set 

H:= R-\K\^ + {ivKf , (11.1a) 
:= 2Di{K^i - ivKg^^) . (11.1b) 

In a field theoretic framework one should provide some further initial data 
for the complete system of equations in the model under consideration; /i and 
u correspond then to the energy and momentum densities of the matter fields 
while Equations (11.1) become constraint equations. Now, we are not assuming 
anything about matter fields, so jj, and v should be thought of as derived from 
the data {M,g,K), as in (11.1). 

Throughout this section we shall be working in a space of spinors which 
carries a representation c of the Clifford algebra associated with the (n + 1)- 
dimensional Lorentzian metric 

ry-- -00^00 + e'®e'= -e^®e^ + g , (11.2) 

with 6^ denoting a co-frame dual to a frame ej. We assume that the c(ej)'s 
are antisymmetric, and that c(eo) is symmetric. The symbol D will be used 
to denote the canonical spinor connection on (M, g) defined by Equation (2.2) 
(and denoted by V there). The connection 

:= Dj + ^K^jc(ei)c(eo) (11.3) 

^Tho results here generalize without any difficulties to spin manifolds of higher dimensions, 
so that the restriction n = 3 is only made for simplicity of presentation. 
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will be called the space- time spin connection on M; V is sometimes referred 
to as the Sen connection. We note that V is not compatible with the positive 
definite metric (•,•)• 

dj{ct>,i,) = (V,-0,V) + (^,V,-V) - (</),K^,-c(e,)c(eo)^) . (11.4) 

However, V is compatible with the Lorentz-invariant (and hence indefinite) 
inner product {4>,tp) := (</>, c(eo)V')- The identity associated with V, which 
replaces the identity (2.6) of Schrodinger-Lichnerowicz, takes the form 

/ (IVVP + i + J^'i^, c(eieo)V)) - I Wl') = / c{ne^)VAi^) , 

JCI JdCl 

(11.5) 

where 

^ := g'^c{ei)Vj . (11.6) 

For sufficiently differentiable {g,Kys, as will be made precise below, the iden- 
tity (11.5) holds in the following circumstances: 

1. If -0 is a iJj^j, spinor field, we may take as O a domain in M with compact 
closure and differentiable boundary dCl; 

2. If is a compactly supported spinor field, then Equation (11.5) holds 
with $7 = M, and dCl = dM; in particular no boundary term is present if 
M has no boundary; 

3. Suppose (M, V) admits a weighted Poincarc inequality, assume dM = 0, 
and let H be the space defined in (8.5). We then have the following: 

Lemma 11.1 Suppose that dM = 0. The function 

Cl(M) 3^^ G(V^) := / (//IV'P + i^*(V',c(e,)c(eo)V)) 

JM 

extends by continuity to a continuous function on M, still denoted by the same 
symbol. 

Proof: For ip,x^ ^ci^) identity (11.9) gives 
G(V)-G(x) = 4/ (-(VV,VV) + (Vx,Vx) + (^V',^V')-(^X,^x)) 

JM 

= 4/ ((Vx- vV',vV' + Vx) + (^x + ^V',^V'-^x) 
< C'llV' + xIIh||V' - xIIh ; 

in the last step we have used the fact that the |c(ei)| are uniformly bounded. 
■ 

Lemma 11.1 implies that the left-hand-side of (11.5) is continuous on IHI for 
any measurable C M. Since compactly supported fields arc by definition 
dense in H, one easily checks, using continuity, that (11.5) holds with Q = 
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M and with vanishing right-hand-side for all E M. (Nonempty compact 
boundaries dM will be considered shortly.) 

We shall say that Mgxt C M is an asymptotically flat end if Mext is diffeo- 
morphic to \ B(0, R) for some R, with 

r~^{9ij - Sij) , dkQij , Kij G L'^{Mey,t) , (11.7a) 
dkgij = o{r-^/^) , K,, = o(r-3/2) , (ll.7b) 
fi, L\M,,t) , (11.7c) 

compare Definition 9.4. Those conditions guarantee that the ADM four-momentum 
of the data set is finite and well defined, as follows from what is said in [21] 
(compare [8]): (11.7a) and (11.7c) guarantee convergence of the mass and mo- 
mentum integrals, while (11.7b) guarantees geometric invariance. 

One of the ingredients of Witten-typc proofs of positive energy theorems is 
the introduction of appropriate boundary conditions on the spinor field in the 
asymptotic regions. In the asymptotically flat case this is straightforward: one 
chooses a ^-orthonormal triad such that, in the coordinate system of (11.7), 

dx''{ei - di) ^r^oo , di (dx''(ej)^ G L^{M^^t) ; 

this is easily achieved by a Gram-Schmidt orthonormalisation of the frame {di}. 
Then a spin frame on Mext is introduced, such that the c(ej)'s are represented 
by constant matrices, as in Section 10. The boundary condition then is that the 
spinor field which will be required to solve the generalized Dirac equation, 

= g'^c{ei)Vj^ = , (11.8) 

asymptotes, as r tends to infinity, to a spinor V'oo which, for r > i? for some 
R, has constant entries in the spin frame above. It is convenient to choose V'oo 
so that ipoo is smooth, and supported in Mext- The procedure is somewhat 
more delicate in the asymptotically hyperboloidal setting; an elegant geometric 
framework for such constructions has been provided in [4, 14] . 

Consider the identity (11.5) with a spinor field = V'oo + X) with x differen- 
tiablc and compactly supported, while dQ = Sr, a coordinate sphere of radius 
R in the exterior region, with R large enough that x vanishes there. A classical 
calculation along the lines of [8] shows that the boundary term in (11.5) is then 
proportional to 

4vrpa('i/'oo ,c(e"e°)V'oo) + o(l) , 

where is the ADM four-momentum of (Afext,^), with o(l) — >■ as i? tends 
to infinity. Passing to this limit we thus have 

/ (I VVI' + 3 (a^IV'I' + ^^'(V', c(eOc(eo)^)) - I^Vl') = 47rp„(V'oo , c(e«e°)Voo) , 

(11.9) 

still for compactly supported x's. But the left-hand-side of (11.9) is contin- 
uous on H, which is shown by a calculation similar to that in Lemma 11.1: Let 
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F{x) denote the left-hand-side of Equation (11.9) with tp = '^^o + X there, let 
Xi G H converge in H to x € so we have 

^(x)-^(xi) = IIxIIh- llxillii 

+2 / (V'Voo ,Vfc(x-X»)) 

JM 

-2 / {&i^oo,nx-xi)) 

JM 

+\ J^i'^'^ ' {l^ + ^''c(efc)c(eo)) (x - Xi)) ■ 

Because VV'oo ^ the first three terms above converge to zero as i — oo. 

The convergence of the final term can be justified by applying the Cauchy- 
Schwarz inequality |(u, Qv)! < -\/ {u, Qu) \/ {v, Qv) whenever Q is positive: 

/ (■0OO , (y" + c{ej)c{eo)) (x - Xi)) 

JM 

< (^Jj^^^°° ' + ^^c(ej)c(eo)) V'oo) j 

X (^^(^ - ^')' + i^^c(ej)c(eo)) (x - Xi))) 

/ \ 1/2 

< C'(II/"IIli + lkllLi)(^llx + XillH||x-Xil|Hj ; 

in the last step of the calculation of Lemma 11.1 has been used. Now, F{xi) = 
F{0), and density implies that (11.9) remains true for any i/j of the form V'oo+X) 
with X G H. 

We are ready now to prove the following version of the positive energy the- 
orem, the regularity conditions of which have been chosen as a compromise be- 
tween those needed for solvability of the Dirac equation (11.8) (c/. Remark 6.1, 
p. 32), those needed for a well defined notion of ADM mass, and those needed 
for a Banach manifold structure for the set of solutions of the general relativistic 
vacuum constraint equations. 

Theorem 11.2 Let (M,g,K) be initial data for the Einstein equations with 
g G W^^^, K G W^ioc' '^^^^ {^i9) complete (without boundary). Suppose that 
M contains an asymptotically flat end and let Pa = ('ti,^) be the associated 
ADM four-momentum}^ If 

fl>\u\g, (11.10) 

then 

m>\p\s, (11.11) 

^"There is a signature-dependent ambiguity in the relationship between p", po and the mass 
m: in the space-time signature ( — , -|-, -|-, -|-) used in this paper this sign is determined by the 
fact that Po, obtained by Hamiltonian methods, is usually positive in Lagrangean theories on 
Minkowski space-time such as the Maxwell theory, while the mass m is a quantity which is 
expected to be positive. 
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with equality if and only if m vanishes. Further, in that last case there exists a 
non-trivial covariantly constant (with respect to the space-time spin connection) 
spinor field on M. 

Remark 11.3 Under the supplementary assumption of smoothness of g and 
it has been shown in [10] that the existenee of a covariantly constant spinor 
implies that the initial data can be isometrically embedded into Minkowski 
space-time, cf. also [66] . We expect this result to remain true under the current 
hypotheses, but we have not attempted to prove this. 

Proof: Suppose that for all ipoo we can establish existence of x = x[V'oo] £ H 
such that ■0OO + X satisfies the Dirac equation (11.8). Equation (11.9) would 
then show that the quadratic form 

ipoo 47rpa(V'oo , c(e"e°)V'oo) 

is non-negative, and the Theorem follows by a standard calculation. The ex- 
istence of X will be a consequence of Theorem 8.8, provided that the relevant 
hypotheses are met. We have 

^ = c{e')Di - ^tTgKc{eo) =^ # = ^ . 

Equation (11.9) with ■0oo = shows that the Schrodinger-Lichnerowicz estimate 
of Definition 8.4, with y = and p = 0, holds both for jO := ^ and its 
formal adjoint = ^. Next, we note that the symmetric part of the 
connection (11.3) is 

= ^K'jc{ei)c{eo) (g) dx^ 

which does not vanish for non-zero X's, but satisfies nevertheless the fall-off 
condition (9.8) by (11.7b). It follows from Theorem 9.5 that the weighted 
Poincare inequality holds. The regularity conditions on the metric imply that 
the requirements of Hypothesis 8.1 with g = g are met: for trivial bundles, or for 
smooth initial data, this is a straightforward calculation, compare Remark 6.1; 
for non-trivial bundles Proposition A. 9, Appendix A, has to be invoked. The 
map (3 of Equation (8.9) is zero, as is the curvature term p in the Schrodinger- 
Lichnerowicz inequality (8.7), by the energy condition ^ > \v\g. From what has 
been said it follows that spinor fields in EI which are also in the kernel of ^ are 
covariantly constant; they then have constant length, and are not in if they 
are non-zero. Theorem 8.8 now shows that for any if)^ there exists a solution 
X G H of the equation 

and the existence of the desired tp follows. ■ 

Let us now turn our attention to manifolds with boundary. We shall say 
that a boundary dM is future-trapped if 

e+:=H+ K{eA, ca) < 0. (11.12) 

A=2,3 



64 



Here H is the mean curvature of dM with respect to an inner-pointing normal, 
while the ca's form an ON basis for TdM. A future-trapped boundary in the 
sense above is future-trapped in the usual sense [36] for a surface in space-time. 
The following result generalises one by Herzlich^^ [38]: 

Theorem 11.4 Under the remaining hypotheses of Theorem 11.2, suppose in- 
stead that M has a diff'erentiable, compact, future-trapped boundary dM. Then 
the conclusions of Theorem 11.2 hold. 

Remark 11.5 One expects that the equality case cannot occur in (11.11), and 
a possible argument could proceed as follows: First, the existence of a covari- 
antly constant spinor implies existence of a non-spacelike, covariantly constant. 
Killing vector field in the associated space-time. Further, if the metric is 
and K is C^, then the space-time metric fulfills the Einstein equations with 
a null fluid as a source [10, Appendix B]; this conclusion is expected to hold 
under the weaker differentiability conditions considered here. By reduction of 
the field equations, this should imply smoothness of the metric. (Alternatively, 
one could assume at the outset that g is and K \s C^, in which case the 
argument presented in the current remark settles the issue). Topological cen- 
sorship results [28] applied to the Killing development [10] of the initial data 
show that the boundary is then the union of a finite number of spheres. Argu- 
ing as in the proof of Theorem 4.6 of [25], the restriction to the boundary of 
the covariantly constant spinor would be harmonic, which is impossible by the 
Hijazi-Bar inequality (2.14) [7,39]. 

Remark 11.6 Past-trapped boundaries are defined by changing the sign of K 
in (11.12); as the remaining hypotheses of Theorem 11.4 are invariant under 
this change of sign, an identical result holds for compact past-trapped dM's. 

Proof: The proof follows closely that of Theorem 11.2, the main difference 
being the need to impose suitable boundary conditions. Indeed, when dM is 
non-empty Equation (11.9) becomes 

/ (I VV'I' + i + i^'i^, c(eOc(eo)V)) - I^V'l') 

Jm 

= 47rpa(V'oo , c(e°e°)^oo) + f (V', c{n)c{e^)V ai/j) 

JdM 

= 4:7rpa{i^oo ,c(e"e°)^oo) + 

/ ^dM^ + k(H- yZiKAAc{n) - KAic{eA)) c(eo)V') , (11-13) 
JdM A 

for, say, continuously differentiable tp's of the form ip = tpoo + with x com- 
pactly supported, and ■i/'oo as in the proof of Theorem 11.2. Further, is 
the Dirac l)()uudarv oixMMtor defined by Equation (2.7), and Equation (2.8) has 

^'^Thc proof in [38] is the rigorous version of an argument proposed in [31]; it also extends 
that argument, as in [31] only marginally trapped boundaries are considered. 
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been used. Finally, is an ON frame on dM with n = ei normal to dM. 
Following [33] we impose the boundary condition 

/C_ := ^(l-e)V' = on^M, (11-14) 

where e := — c(n)c(eo). We then have 

{'(p,c{eA)c{eo)tp) = -{^p,c(eA)cieo)cin)c(eo)^J) 

= -(V',c(eo)c(n)c(eA)c(eo)V') 

= -(c(n)*c(eo)V,c(eA)c(eo)V') 

= (c(n)c(eo)V',c(eA)c(eo)V') = -{'^P,c{eA)c{eo)^p) , 

which shows that the last term in the last line of Equation (11.13) vanishes. 
Then 

and e* = e, so 

= -{'ip,e%M'>P) 

which shows that the first term in the last line of Equation (11.13) vanishes. 
Next, 

A A 

which shows that the sum of the second and third term in the last line of 
Equation (11.13) gives a non-positive contribution when dM is trapped. When 
M ^ W\gi V'oo = 0, and (11.14) holds, from Equation (11.13) we obtain 

/ (|VV'|' + i(/x|^P + i^^(^,c(e,)c(eo)V)))-i/ 0+(V,V')=/ I^V'P , 

JM JdM Jm 

so conditions (11.10), (11.12) give 

/ IVV'I^ < / i^vP , (11-15) 
Jm Jm 

for all V e C'i(M) which satisfy (11.14). Define 

M,C- ■■= G H, /C-V' = on dM} , (11.16) 
where H is defined in (8.5). If we let 

G'W:= [ {l{fi\^Pf + ly\^P,c{eMeo)^P)))-k{ 0+{^,i^) , (11-17) 
Jm JdM 

then the calculation of the proof of Lemma 11.1 shows that G'{'ip) can be ex- 
tended by continuity to a continuous function on Hx;_ . The boundary integral 
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in (11.17) is continuous, so the volume integral is also continuous on M.jc_ . which 
implies that Equations (11.15) holds for all -0 € ]HIa:_. Lemma 8.6 establishes 
the Schrodinger-Lichnerowicz estimate (8.7) with p = for (^, /C_). 

As explained in Section 10, the boundary value problem determined by 
Equation (11.14) belongs to the family of problems considered in Theorem 8.8. 
If both M and dM are simultaneously parallelizable, and if the metric g is 
a product near dM, then the regularity conditions of Theorem 8.8 are met 
by hypothesis; the general case is handled by Propositions A. 9 and A. 10, Ap- 
pendix A. Repeating now the arguments of the proof of Theorem 11.2 gives 
the non-negativity of po- ■ 

It is expected that the positivity statement of Theorem 11.4 can be strength- 
ened to the so-called Penrose inequality when trapped boundaries occur. This 
question remains wide open, except in the special^^ case Kij = [15,42]. An 
interesting related inequality has been, essentially, proved by Herzlich [37]; how- 
ever, the arguments of that last reference do not include a sufficient justification 
of existence of the required spinor field, except in the rather special case of a 
smooth metric which is a product near the boundary, as analyzed by Bunke [16]. 
Here we fill this gap and establish the following: 

Theorem 11.7 Let {M,g) be a complete Riemannian manifold with g G W^^^, 
and suppose that M has a boundary dM diffeomorphic to S"^, with non-positive 

inwards pointing mean curvature. Suppose that the curvature scalar R{g) of the 
metric g is non-negative, and that M contains an asymptotically fiat end with 
mass m. If a is the dimensionless quantity defined as 



Moreover, if the metric is smooth, then equality is achieved if and only if 
{M, g) can he isometrically embedded in the Schwarzschild space-time with mass 
^/Arm^^M)Jm^. 

Remark 11.8 If M is the union of a compact set with a finite number of asymp- 
totically fiat ends, then a > 0. 

Proof: The details of the argument follow closely those of the proof of The- 
orem 11.4, the pointwise boundary conditions (11.14) being replaced by the 
spectral boundary conditions (2.10) with /C given by (2.13); compare the dis- 
cussion of Section 2, as well as that in the paragraph following Equation (10.7). 
The main elements missing in the arguments of [37] are provided by the bound- 
ary regularity results of Section 5; those arc the key to the proof of Theorem 8.8. 

^^Similarly to Theorem 11.7, for the results in [15,42] it actually suffices that ii > <4> /j > 
— l-f^lg + (tTgK)'^, and that tr/j/sT vanishes on dM, where tr^ is the trace of the restriction of 
K to dM. 




inf 

/ec-(M),/^o 




then 




67 



The reader is referred to [38] and [23, p. 679] for the analysis of the equality 
case, cf. also [9]. ■ 

Following [31,34], let us pass now to inequalities with an electric charge 
contribution. A set (M, g, K, E, B) will be called initial data for the Einstein- 
Maxwell equations if (M, g) is a three dimensional Riemannian manifold, K is 
a symmetric tensor on M, while E and B are vector fields on M; as before, 
this is a serious abuse of terminology, as we are not requiring any constraint 
equations to be satisfied. Given such a triple we set 

div E := DiE' , 
div B := DiB' , 
^^■.= R- \K\l + {tvKf - 2\E\l - 2\B\l , 
Ui := 2Dj{Kh - trKSj) + ieijkE^B'' . (11.18) 

Here D is the Levi-Civita connection associated with the metric g. In a general 
relativistic context, div E is the electric charge density, div B is the magnetic 
charge density (usually zero, whether electro- vacuum or not), ^ is the energy 
density remaining after subtracting the electro-magnetic contribution, and is 
the left-over matter current; /x, div E, div B and Jj vanish when the Einstein- 
Maxwell constraint equations hold. 

Let a new connection V be defined as 

Vi := A + ^i^iic(e^>(eo)-^c(£;)c(ei)c(eo)-^ejfc^S^c(e*^)c(e^)c(eO ; (11.19) 

The connection V will be called the space-time Einstein-Maxwell spin connec- 
tion on M. V is again not metric compatible, with symmetric part given 

by 

= (^Kijc{e^)c{eo) - Eic{eo) - ^e,fc^5^c(e'=)c(eOc(eO) ® e' , (11.20) 

where, as before, 6^ is the co- frame dual to ej. In this context the asymptotic 
flatness conditions have to be complemented by conditions on E and B: we 
shall require 

E,BeL\M,,i), divS,divSGLHMext) , ^ = ©(r"^) , B = o{r-^) . 

(11.21) 

Following an argument proposed by Gibbons and Hull [34] we have: 

Theorem 11.9 Let {M,g,K,E,B) be initial data for Einstein-Maxwell equa- 
tions with g G M^ioc' K,E,B G W^^^ , with {M,g) complete (without boundary). 
Suppose that M contains an asymptotically flat end M^^t with E and B satis- 
fying the fall-off conditions (11.21) there, and let Pa = {m,p) be the associated 
ADM four-momentum. Let Q and P be the total electric and magnetic charge 
ofM^^, 

Q = lim -i- / E'dSi , P = hm / B'dSi . 

i?— >oo 47r Jr=R i?— >oo AtT Jr=R 
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If 

> ^JW\g + |div + |div B\^ , (11.22) 

then 

m>^\p\l + Q'' + P^ , (11.23) 

where \p\s ^ V^ip^; ^^^^^ equality if and only if there exists a spinor field on 
M which is covariantly constant with respect to the Einstein-Maxwell space-time 
spin connection (11.19). 

Remark 11.10 Under the hypothesis of smoothness of the metric, Tod [61] has 
found the local form of the metrics which admit covariantly constant spinors 
as above; however, no classification of globally regular such space-times is 
known. It is expected that the only singularity-free solutions here have vanish- 
ing Maxwell field, or belong to the standard Majumdar-Papapetrou family (c/., 
e.g., [27]). It would be of interest to fill this gap. 

Remark 11.11 Charged matter might violate (11.22); however, there might 
exist a constant a G (0, 1) such that 

H > + a2(|div E\'^ + jdiv B|2) . (11.24) 

Replacing in (11.19) the fields E and B by aE and aB, an essentially identical 
argument leads to 

m> ^[pl2 + a2(Q2 + p2) . (11.25) 
Proof: For the connection (11.19) the identity (11.9) becomes [34] 

JM ^ 

+\{'^, + z/*c(ei)c(eo) - div Ec{eQ) — div Bc{n)c{e2)c{e3)) 

= 47r(V'oo , [p„c(e°)c(eO) + Qc(eo) - Pc{n)c{e2)c{e3)] ^^o) , (11-26) 

again for of the form i/joo + X) with compactly supported x's- The Dirac 
operator 

is not formally self-adjoint, we have instead 

# = c{e')Di -\(tVgK - c(£;))c(eo) + ^ejuB^ c{e'')c{e^) . 

This shows that the adjoint of ^ coincides with ^ modulo the replacement 

B~^ -B . 

The arguments follow now the previous ones, basing on the identity (11.26). We 
simply note that (11.22) implies non-negativity of the quadratic form appearing 
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in the second line of (11.26). Similarly, positivity of the quadratic form defined 
by the third line of Equation (11.26) implies (11.23). Some comments arc in 
order here, related to the fact that L is not formally self-adjoint when the 
magnetic field does not vanish. Since we are not assuming interior compactness 
of M, M could have other asymptotic regions in which B could grow in an 
uncontrollable way, so that will not map EI into I? . Now, C) differs from L 
by a change of the sign of iJ, which implies that L'^ also satisfies a Schrodinger- 
Lichnerowicz identity with a connection in which B is replaced by — B. The 
arguments already given show that the weak equation £^0 = has no I? 
solutions, and Corollary 8.8 provides the desired isomorphism property of £. ■ 

In the presence of boundaries we have: 

Theorem 11.12 Under the remaining hypotheses of Theorem 11.9, suppose in- 
stead that M has a compact future-trapped boundary dM. Then the conclusions 
of Theorem 11.9 hold. 

Proof: This is a repetition of the argument of the proof of Theorem 11.4; one 

imposes again the boundary condition (11.14), and we only need to check that 
(11.15) still holds. This is indeed the case, which is established as follows: the 
electromagnetic field leads to a supplementary contribution 

/ (V', [E'c{eo) - B'c{e^)c{e^)c{e^)\ ^)ni 
JdM 

to the boundary integral (11.13). When (11.14) holds we have 
(V',c(eo)'0) = -(V',c(eo)c(n)c(eo)'0) 

= -{c{n)i>,^p) 

= ~{c{n)c{eo)c{eo)i>,^p) 

= {c{eo)c{n)c{eo)ip,ip) 

= -{c{eo)ip,ip) = -(V',c(eo)V') , 

hence 

(^,c(eo)V) =0. 
Similar manipulations show that Equation (11.14) implies 

(V',c(ei)c(e2)c(e=^)V')=0, 

and the result follows. ■ 

We finish this section by noting that positive energy results follow by iden- 
tical arguments for asymptotically hyperboloidal manifolds [4, 25, 26, 29, 33, 52, 
62,67]; here Theorem 9.10 should be used instead of Theorem 9.5. The defini- 
tion of mass in that case is considerably more delicate, we refer the reader to 
[25, 26] for details. 
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A Fields on manifolds of W^'^^'^ differentiability class 



Consider a smooth manifold M; on such a manifold one can define in a geomet- 
rically invariant way tensor fields which are of C°° differentiability class, or of 

class, or of class. For example, one says that a tensor field is of Wl^f 

class if there exists a covering of M by coordinate patches such that the coordi- 
nate components of the tensor in question are in in each of the coordinate 
patches. Since the transition functions when going from one coordinate system 
to another arc smooth, this property will be true in any coordinate system. 

Let, now, (M, g) be a smooth manifold with a pseudo-Riemannian metric g 
which is of differentiability class. For various arguments it is convenient to 
use local coordinate systems which are adapted to the metric, such as geodesic 
coordinates, or harmonic coordinates. In this case the transition functions to 
the adapted coordinate system will not belong to the original smooth atlas on 
M in general. At this point there are two strategies possible: either to enlarge 
the atlas on M to contain those new coordinate systems, or to ignore this issue 
and try to analyze the problems that arise on an ad hoc basis. For nearly 
all of this paper the ad hoc approach, working entirely within a structure 
on M, is quite adequate. However, the proofs of Theorems 11.7 and 11.12 
require the existence of approximately Gaussian coordinates near a boundary 
or near a hypersurface of M {cf. Lemma A. 10 below). Direct construction 
of such coordinates with respect to the metric g produces a coordinate 

change which is not C°° , which forces us to analyse the problems involved when 
constructing systematically manifolds of Wy^~^^'^ differentiability class. For this 
reason we will present such a construction here. For technical reasons we shall 
always assume that 

p£[l,oo], keN, kp>n; (A.l) 

these restrictions are more than sufficient for our purposes. Generalising the 
condition k e N to k e M"*" would require an analogue of Lemma A. 2 for 
non-integral k,£, which seems not to be available. Condition (A.l) and the 
Sobolev embedding W^^'^{^) C C^{i}) (for appropriately regular open do- 
mains C M") mean that we will consider only manifolds which are at least of 
differentiability class. 

Consider, thus, a connected paracompact Hausdorff manifold M of dif- 
ferentiability class. We shall say that M is of W^^^^'^ differentiability class if M 
has an atlas for which all the transition functions arc of W^^^'^ differentiability 
class. Unless indicated otherwise, the Lebesgue measure in local coordinates is 
used. 

A tensor field with components which are C°° with respect to some coordi- 
nate chart (belonging to a C°° sub-atlas of the Wj^^^'^ atlas), will not generally 
have smooth components in all W^^^^'^ charts. In this situation a tensor 
field "comes equipped" with a preferred atlas of coordinate charts in which it 
has smooth coordinate components. This is a priori the case for tensor fields 
of any differentiability class on W^^^'^ manifolds, and it is of interest to single 
out those classes of tensor fields, the coordinate components of which will be 
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of a prescribed differentiability class in every coordinate system of the W^^'^ 
atlas on M. Differentiability classes of this type will be referred to as invari- 
antly defined. Our next result describes some such classes of tensor fields. It 
is convenient to introduce the following notation: let x,y eM., we shall write 
X >* y if the following holds: 

x>^y^ { ./.n' A.2 

lx>y,ny<0. 

(We note that for x > the only value of x at which does not coincide 
with ">" is X = 0.) In this notation the the Sobolev embedding theorem can 
be stated as: 



Proposition A.l Let {M,g) be a W^^'^^'^ manifold, kp> n, pE [l,oo]. 
1. Let {£, q) be such that the Sobolev embedding 



loc loc 

holds (which is equivalent to the condition 

->*- + , A.4 

q p n 

with >* defined in (A.2)). Then the space of W(^^ scalar fields on M is 
invariantly defined. 

2. Let {£, q) be such that the Sobolev embedding 

holds (which is equivalent to the condition 

- >* - + .) A.5 

q p n 

Then the space of W^^^ tensor fields on M is invariantly defined. 

Proof: Point 1 is a straightforward consequence of the following Lemma: 

Lemma A.2 Let Q,,U C and let i/j : Q ^ U be a diffeomorphism such 
that tjj G l¥j^^^'^(r2;M'*), kp > n. If {£,q) is such that the Sobolev embedding 
Wi^^^'^ C Wl^^ holds, cf. Equations (A.4) and (A.2), then for all F e w/„'^(W) 
we have 

Remark A. 3 In [13, 58] some partial results can be found concerning sharpness 
of this result. 
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Proof: We have, for < |a| < £{< k + 1), 

d''{Fo,p) = ^C{au . . . ,am)d''^tp ■ ■ ■ d''"^tpF^"'K i; , 

where the sum is taken over sets (cti, . . . , a„i) satisfying ai + • • • + am = ol, 
with \ai\ > 1. For any compact K (Z^ it foUows that 

1 111 

+ • • • + 7 + 7 < - ■ (A.7) 



s{ai) s("m) * 

Here we have used the generahzed Holder inequahty, and the change of vari- 
ables theorem to pass from o V'||Lt(A:) to H-^*-™'' By Sobolev's 
embedding we have F^"^^ G L'^{ip{K)) for all r satisfying 

I I m- e 

- >* - + . A.8 

r q n 

Consider, first, those terms in (A. 6) for which the right hand side of (A.8) is 
positive (if any). Let r be defined by Equation (A.8) with replaced by =. 
Set 

' " ' ; (A.9) 



since kp > n we have 



s{ai) n 



1 1 k \ai\ — 1 , , 

> + — , (A.IO) 



s{ai) p n n 



and Sobolev's embedding theorem implies that ||9'"*''i/'||2;^s(ai)^;^') is finite. With 
this choice of r and of the s(ai)'s we have 



y- 



1 \a\ — m 1 m — £ 
+ - = — + - + 



s{ai) r n q n 

\a\-e 1 1 

+ -<- , 



n q q 

so that those terms will give a finite contribution to the right hand side of (A. 6) 
by setting t = r. 

Consider, next, those terms in (A. 6) for which the right hand side of (A.8) 
vanishes. If all the aj's have length one the term in question will give a finite 
contribution to the right hand side of (A. 6) by setting t = q. If one of the a^'s, 
say CKi, has length large than 1, for z > 2 we choose the a^'s as in (A.9), while 
we set l/s(ai) = (|ai| — 1)/?^ — e > 0, with < e so chosen that (A.IO) still 
holds, e < {kp — n)/2pn. Choosing 1/t = e will lead to a finite contribution in 
(A.6). 

It remains to consider those terms in (A.6) for which the right hand side of 
(A.8) is negative. In this case we set t = oo, and 

N v^-M<- fAin 



s{ai) £q ^ s{ai) £q q 
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By Sobolev's embedding ||5°'V'llL''(«2)(i<') '^^^^ t)e finite when (A. 10) holds. Now 
Equation (A. 10) with s{ai) defined by (A.ll) is equivalent to 

t \q nJ p n 

The right hand side of (A.12) is negative. If the left hand side is positive or 
vanishes there is nothing to check. If both sides are negative the worst case is 
obtained with \a.i\ = i, and the inequality holds when (A. 4) is an inequality. 
The simple analysis of the case of equality in (A. 4) is left to the reader. ■ 

Before returning to the proof of Proposition A.l we need one more Lemma: 

Lemma A. 4 Let < m < I < k, q,p £ [l,oo], kp > n. Suppose that {i,q) is 
such that the Sobolev embedding W^^^ C W^^^ holds, cf. Equations (A. 5) and 
(A. 2). Then the product map 

Wtr'^xwtl 3 if,g) -^fgeWl--''^ 

is continuous. 

Proof: For any < \a\ < £ — m the Leibniz rule gives 

5"(/5)= E C^i,a2d"'fd°''g , 

SO that on any compact set if C M the Holder inequahty gives 

P"(/5)IIl.(x)<C Y1 , (A.13) 

ai+a2=oc 

^ + -7^ < - ■ (A.14) 



s(ai) s(a2) ~ q 
Let 

/ s 1 \ai\ — k + m 

a(ai) = - H , 

p n 

1 \a2\ -i 

a{o(2) = - H • 

q n 

By Sobolev's embedding we will have 9"^/ G L^("i)(E:), d^'^'g G U^°'^\K) when 

1 



>* a{ai) . (A.15) 



s{ai 

We have the following cases: 

• If a{Q.i) < and 0(02) < we set s{ai) = q/2, and we obtain 

II'9"VIIl»{c«i)(k)II^"^5'IIl''(«2)(k) ^ C\\f\\w''-^-p{K)\\9\\w^.i{K) (A.16) 
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• If a{ai) > and 0(02) > we set l/s{ai) = a{ai) so that (A. 15) holds, 
and we obtain 

1 1 1 1 \a\-k-e + m 1 1 k 1 

+ r = - + - + — <- + <-, 



s(ai) 5(02) q P n ~ q p n q 

since kp > n, so that (A. 14) holds. We note that (A. 16) is again satisfied. 
If a{ai) = and 0(02) > we have 

\a2\ |a| — |ai| \a\ — k + m 1 



n n n p 

so that 

1 1 \a\ — k — i + m 1 1 k 1 
a 02) = - + - + ^ <- + < - . A.17 

q p n q p n q 

Let e be any number satisfying < e < {kp — n)/np, set 

1 e 1 , ^ e 

a{a2) + 



s(ai) 2 ' s(a2) 2 ' 

Decreasing e if necessary we will have 5(02) > 0. Then (A. 14), (A. 15) 
and (A. 16) hold by the calculation in Equation (A.17). A similar analysis 

takes care of the case a(ai) > and 0,(02) = 0. 

• If a(ai) < and 0(02) > we set s(ai) = 00 and 5(02) = q- 

• If 0(0:1) > and 0(02) < we set s{ai) = q and 5(0:2) = 00; (A. 14) 
obviously holds, while 

1 ,,11 A;-m-|oi| 1 1 k-i 
— 0(01) = 1 > h 



s(oi) q P n ~ q P n 

which is non-negative by (A. 5), hence (A. 15) and (A. 16) hold again. 

This establishes that fg G M^ioc'"'^- "^^^ continuity follows immediately from 
the inequality 

\\fg\\w'^-m,q(^K) < C!\\f\\iYk-m,p(_K)\\g\\wi'i{K) 

which has been established during the proof. ■ 

We can pass now to the proof of point 2 of Proposition A.l. Let two coor- 
dinate systems on M be given related to each other by a map ip € ^itc^^'^' 
X = V'"^- Let i"^'""''/3i.../33 and i'^^'"'^''vi..Ma be the coordinate components of a 
tensor field t, with G W^''^. We have the transformation rule 

t^-^^.,....M = t--"^/3.../3.(^(x))^(x) . . . ■ ■ ■ ^(^(^ 

(A.18) 



75 



Now the matrix o -0 is the inverse matrix to^^, so the components of the 
former are rational function of those of the latter. We recall the Gagliardo— 
Moser-Nirenberg inequalities {cf., e.g. [41, Corollaries 6.4.4 and 6.4.5]) 

V f,g e w^'^nL- WfgWwe. < Ci(||/|Uoo||5||^,„ + (A.19) 

V / G W'''^ n \\F{f)\\we. < C2(||/||loc)(1 + 11/11^...) (A.20) 

for some / and g independent constant Ci, and for some constant C2(||/||l°°) 
depending upon / only through its norm. Here one assumes that F is a 
smooth function of its argument which is allowed to take values in R^, and the 
integrals are taken over compact sets. Equations (A.19)-(A.20) show that the 



V''s are W^'^ functions of their arguments. Lemma A.4 implies that 



It follows that the right hand side of Equation (A. 18) is of the form 

toi;A, te Wt , ^ e , A e W^f . 

Lemma (A. 2) implies that toip £ ^loc' Lemma A.4 with m = shows that 
the right hand side of Equation (A. 18) is in W^^'^, as desired. ■ 

We wish to extend the above discussion to spinor fields; this requires the 
introduction of orthonormal frames, and hence of the metric. Consider, then, 
a VFj^^^'^ manifold M with a strictly positive definite symmetric two-covariant 
tensor field g. We shall say that (M, g) is a pseudo-Riemannian VFj^^^'^ mani- 
fold if M is a VFjq^^''^ manifold and if 5 is a pseudo-Riemannian metric of W^f 
differentiability class. This is an invariantly defined notion by Proposition A.l. 
Before proceeding further we note the following: 

Proposition A. 5 Let (M, g) be a Wl'J^^'^ manifold with a pseudo-Riemannian 
metric of W^^^ differentiability class, kp > n, p £ [1, 00]. Then the following 
hold 

1. In any coordinate system in the W^^^'^ atlas the Christoffel coefficients 
r^fc satisfy 

k-l,p 



2. The Riemann tensor is of W^^'^'^ differentiability class. 

3. The curvature scalar R = g^-'R^ikj is ofW^~^'^ differentiability class. 

4- Assume that i > 1 and suppose that {i, q) is such that the Sobolev embed- 
ding W^^f C W^^^ holds, cf. Equations (A. 5) and (A. 2). Let t be a tensor 
field of VFio'p differentiability class, then for any vector field X G W^^^ we 
have 
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Proof: 1. By definition of the T^j.'^ h.aYe 



pi a (^9ei dgik _ dgjk \ 

By (A. 19) we have g^^ G W^'^, by definition the derivatives of the metric are in 
W^~^'^, and lemma A.4 with m = and {£, q) = {k — l,p) gives the result. 
2. By definition of the curvature tensor we have 



T>i J'' I pi pm pi pm 

~ 'dx'^ "dx^ '^^ ""^ ' 

Lemma A.4 with m = 1 and g) = {k — l,p) shows that the product terms 
above are in W^jo"^'^, and the result follows from point 1. 

3. This claim follows immediately from point 2, from g^^ G VFj^jf and from 
Lemma A.4 with m = and {i,q) = {k — 2,p). 

4. Let i be a tensor field of W-^^^ differentiability class, in local coordinates 
we have 

The r's are in W^^~^'^ by point 1, thus the product terms are in W(^^''^ by 
Lemma A.4 with m = 1. The claim about X^'Vit follows again from Lemma 
A.4. ■ 

Let {M,g) be a Wjo^^'*' pseudo-Riemannian manifold and let OAfM,g be 
the bundle of ^(-orthonormal frames on M. We can equip OJ\fM,g with a Wj^jf 
structure by considering only those ^(-orthonormal sets of vector fields which are 
all of W^J^ differentiability class. Let us start by showing that the set of such 
(locally defined) frames is not empty. On O, the domain of a coordinate system 
(x*), we can construct a (/-orthonormal frame Cj = ej^djdx^ by performing a 
Gram-Schmidt orthonormalisation of the basis By construction the 

coordinate coefficients e^* of the vector fields Cj are smooth functions of gij (at 
least on a neighborhood of the range of values taken by gij), where the g'ij's 
are the coordinate coefficients of the metric g, g = gijdx'^dxK Since kp > n, 
(A.19)-(A.20) applied to the e^* considered as functions of the gij shows that 
the vector fields ej are indeed of W-^J^ differentiability class, as desired. 

The following shows that the W^Jj^ structure of OMu^g is an invariantly 
defined property of a W^^'"^ pseudo-Riemannian manifold: 

Proposition A. 6 Any two (globally or locally defined) g-orthonormal frames 
of W^^f differentiability class are related to each other by a 0{n) -rotation of 
^\oc diff^f^ntiability class. 

Proof: Consider two locally defined ^r-orthonormal frames Cj and fi, i = 
l,...,n = dimM, of W^^f differentiability class. In particular each of the 

gj and /j' is a vector field of wl'^f differentiability class, which is invariantly 
defined by Proposition A.l, so that it is sufficient to prove the result in any 
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coordinate system in the W^^^'^ atlas on M. In such a coordinate system {x*} 
we can write ej = e/^, fj = //^i for some functions e/,// G ^lof- Since 
both frames are orthonormal there exists an 0(n)-valued function Wi^ such 
that 

ei = wifj . (A.21) 

It follows that 

where f^^. is the matrix inverse to fj^. We have f^/^ G W^^f by (A. 20), thus 
Wi^ G Wj^f by (A.19), hence the result. ■ 

Now suppose that AI has a spin structure, namely a Spin-principal bun- 
dle ^ which double-covers the principal bundle ^ of g'-orthonormal frames of 
{M,gy. 

0^Z2^^^^^M. (A.22) 

We note that the obstruction to the existence of such structures is purely topo- 
logical, (c/., e.g., [45, Chapter II]) and therefore independent of the choice of 
the metric and differentiability class. A bundle of spinors Y = ^ Xt V is a 
vector bundle associated to ^ and a representation 

T : Spin ^ End(F) , 

for some finite-dimensional vector space V. A choice of ^(-orthonormal frame 
e = (cj) of H^iq'^ differentiability class defined on an open set O C M determines 
a local section of This lifts to a section of which in turn is associated 
with a local orthonormal frame 4> = {4>i) in Let U be another open set 
with a g orthonormal frame e' = (e'-), so by Proposition A. 6 there exists an 
0(n) valued map w = {w,p) of VFj^jf differentiability class such that the frames 
e, e' are related by Cj = Wi^e'^ on O n The map w : O CiU ^ SO lifts to 
w : O nU ^ Spin. This lift is not unique, but the possible lifts differ only by 
a fixed nontrivial element z of the centre Z(Spin) ~ Z2. The corresponding 
spin frames (j),(f)' are related by = T{w)(j)' or T{zw)(j)' . Analyticity of the 
local inverse vr^-'^(-) and the inequality (A. 20) show that T(vj),T(zw) G W-^^f, 
and it follows from (A.19) that the spin frames (j) on U and cp' on O are VFj^jf 
compatible onU DO. This establishes the following result: 

Proposition A.7 Let (M,g) be a VFj^^^'-'' spin manifold with a pseudo-Riemannian 
metric of Wj^'^ differentiability class. Then every spinor bundle carries a nat- 
ural W^J^ differentiable structure. 

An argument similar to that of Proposition A.l shows: 

Proposition A. 8 Let {M,g) be a W^^^^'^^ spin manifold with a pseudo-Riemannian 
metric ofW^^^ differentiability class, kp> n, p ^ [^,00]. Let (i,q) be such that 
the Sobolev embedding 

holds, cf. Equations (A. 5) and (A. 2). Then the space 0/ Wj^'^ spinor fields is 
invariantly defined. 
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To proceed further, we recall the definition of the covariant derivative of 
a spinor field. Let V, €{, O C M and 7^ be as before, and let (p = 
be the spinor frame corresponding to the orthonormal frame e = (cj). This 
defines a preferred local spin frame, with respect to which the Clifford action 
is represented by locally constant matrices 7j. Let be a spinor field over 0\ 
the spinor covariant derivative of (j) is given in terms of the orthonormal frame 
connection matrix LOij{ek) = g{ei,V e^^j) and the spinor frame components 
ip = by (2.2), and the Dirac operator of V on S is defined by 

= y Ve,V • (A.23) 

Proposition A. 9 Let {M, g) he a W^^^^'^ manifold with a pseudo-Riemannian 
metric of W^^^ differentiability class, kp > n, p ^ [l,oo]- Then the following 
hold 

1. Let be any g -orthonormal frame ofW^^^ differentiability class, then the 
spin connection coefficients defined as Vg^.V' = ^k{'^) +^k'^ satisfy 



2. If {i,q) is as in Proposition A. 8 with i> 1, and if X is a vector field of 
W'^'P differentiability do 
to spinor fields: 



W'^'P differentiability class, then Vx maps continuously VFj^'^ spinor fields 



In particular the Dirac operator maps continuously W^^^ to 

Proof: To prove point 1 choose a spin frame in which the c(e*)'s are point 
independent matrices. Then 

Wfe = -^c{e^)c{e^)ujij{ek) , 

with ojij{ek) = g{ei,V ek^j). The claim that u}ij{ek) G l^itc^'^ follows immedi- 
ately from point 4 of Proposition A. 5 and from Lemma A.4. The result in any 
spin frame follows from the transformation rule of the connection coefficients 
under changes of frames and from Lemma A. 4. The proof of point 2 follows 
that of point 4 of Proposition A. 5 and will be omitted. ■ 

Given a smooth metric in a neighbourhood of a compact boundary dM, 
geodesies normal to the boundary determine a diffeomorphism of F x 7 with 
a neighbourhood of dM ~ Y, such that in adapted coordinates v = {y^,x), 
X G [0, xo), y"^ G Ui we have 

= 1, = 0. (A.24) 

The diffeomorphism determines a tubular neighbourhood of dM and the re- 
sulting coordinates are called Gaussian coordinates. If the metric has only 
low differentiability then uniqueness of the geodesic equation may fail, and the 
existence of Gaussian coordinates becomes problematic. However, for our ap- 
plications it is sufficient for (A.24) to hold only approximately near dM, in 
which case we may rely on the following result. 
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Proposition A. 10 (Almost Gaussian tubular neighbourhood coordinates for 
dM). Lei G N, £ G N U {0}, and suppose {M,g) be a T^j^^^'^ Rieman- 
nian manifold with metric g G Wjq^(M), {k — £)p > n. Let Y C M be a 
compact connected component of the boundary of M with Y of W^'^^'^ differ- 
entiability class. There is a neighbourhood O ofYdM and x G W^^^'^{0), 
and a diffeomorphism O Y x I , I = [0,.xo), which determines coordinates 
= (y^, x) G O such that F n O = {a; = 0} and 

g{dv\ dv^) = g'^ G W''^p{0) , (A.25) 
g{dx ,dx)-l = 0(x^+'^), (A.26) 
g{dx,dy^) = 0{x^+'')., (A.27) 

for some a > 0. 

Remarks: 1. A similar result for C*^'^ metrics, fc > 1, follows from [3, Ap- 
pendix B] . 

2. Similar results hold for pseudo-Ricmannian manifolds provided Y is 
non-characteristic, and when y is a hypersurface in M. 

Proof: If = go we can use Gauss coordinates near Y, and the result follows. 
Suppose thus that A; < oo, let x be any defining function for Y and let Oa be any 
conditionally compact coordinate neighborhood of Y, with g^^ = g{dv^,dv^) G 
W'^'PiOa). Passing to a subset of Oa if necessary without loss of generality we 
may assume Oa « [0,xo) xUa,Ua C Y. Coordinate systems of this form will 
be called cylindrical. 

We construct a suitable M/'^+^'P coordinate change {y^,x) in Oa X 7 by 
noting first that dx = ^dx + -^dy^, dy^ = %^dx + ^jgdy^. Thus if x is 

also a boundary coordinate, so x(y, 0) = and = on y, then the metric 
coefficients satisfy 

g^^ = g{dx, dx) 

XX I 9x I dx dx 
dx J dx dy^ dy^ dy^ 

g)^^^ onY = {x = 0}, 



= g{dx,dy^) 



dx ( dy^ dy^ 



+ ^ ) on y 



dx \ dx dy 

Since (f^ G W^'^, the restriction (f^\Y lies in the Besov space A^'^^^^(y) (see 
[59, §VI.4.4], or [43, Theorem VII. 1]) and there is an extension x = x{y,x) G 
W'^^^'P{Ua X I) satisfying the conditions 

x(y,0)=0, ||(y,0) = (5^-)-V2(y,o), 

for y = (y^) G Oa ([59, §VI.6], [43, Theorem VII.3]). This implies 

= g{dx, dx) G W^'P{Ua x /) 
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and (f^ = 1 when x = 0. Similarly, there is f^{y, x) G W^'^^''P{Ua x I) such 
that 

f\y, 0) = 0, ^(y, 0) = - g^^'/g^^iy, 0) , 
so the coordinates (y, x), y^{y, x) =y^ + f^{y, x) also satisfy 

g^^{y,0)=g{dx,dy^) = 

on y, since dy^ = dy^ — j g^^ dx. 

As shown in Proposition A.l, in the new coordinate system (v*) = {y^,x) 
we still have g{dv^,dv^) G W^'^{Oa), so that by embedding theorems the metric 
coefficients are cr-Holder continuous on Oq, for some cr > 0, and 

- 1 = 0(x<^), g^^ = 0(x<^) . (A.28) 

Let ^, G VF*^"'"^'^(C>q;), and consider the effect of the change 

x = x + (t){y,x), (A.29) 

y'' = / + /^(y,x), (A.30) 
where 0) = 0, /"^(y, 0) = 0. Then 

g^^ = g-- + 25^^^ + ff^^'^l^ , (A.31) 



9^'^ ^ g{dx, dy^) = g-^ + 5-^^ + 5^^^ + g'^^'^- (A.32) 



Suppose that for some ^ > we have 

g^^-l = 0(x^+'^), = 0(x^+'^) . (A.33) 

This holds for £ = by (A.28) and we establish the general case by induc- 
tion. Again by restriction and extension results [43, 59] there exist (p, G 
satisfying 



dx 



1 cfg 



\f. „xx 



{y,o) 2 dx^ 

ge+lfA ^ QigXA 

' ',0) 



0) 



(2/,o) dx^ 



while all the lower order x-derivatives of c/),/^ vanish at x = 0. Passing to 
coordinates {y,x) on a (possibly smaller) cylindrical neighborhood Oa, one 
finds from (A.31)-(A.32) that (A.33) stiU holds and moreover, 

/^-l = 0{x^+^+''), 
= 0(ar^+i+<^), 

Dropping bars one finds that (A.33) holds with £ replaced by £ + 1, and the 
induction step is complete. 

Finally we show that the local charts can be combined to form a tubular 
neighbourhood diffeomorphism. It follows from (A.31) that if x, x both satisfy 
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(A. 28) and vanish on Y, then g{dx,dx) = 1. In particular, by combining the 
functions Xa from each of the local coordinate charts Oq, using a subordinate 
partition of unity (pa, the function x = T,a<pa^a satisfies x G W^^^'P{Y x /) and 
X = 0, g^^ = 1 on y. 

In order to construct a diffeomorphism with y x / we need to construct 
a similar averaging of the coordinate functions. Fix a smooth embedding 
$ : y ^ R-'^ and let n$(y) : M C ^(^) be the orthogonal projection 

in M.^ from a tubular neighbourhood J\f back to ^{Y). Let ya = (ya) '■ C 
Y W^'^ denote both the coordinates of a C°° chart on Y, and their natural 
extension to y^ = (ya) : Oa = Ua >^ I C Y x I ^ M""^ Let y^ = : 
W^~^ be the functions constructed above, so there is a neighbourhood Oa C 
containing Y = Y x {0} such that y^^ o y^ : Oa Y. Choose a finite covering 
Oce oiY X I with subordinate partition of unity 4>a and define :Y x I ^Y, 

^{p) = o n^^Y) {ZaMp)Hya ' ° Up))) ■ 

Since y^iy-,^) = y^, *|yx{o} = Id and defines a diffeomorphism of 

y X /. Now for any C°° chart y = (y^) on Y, y^ := y^ o defines a 
chart on y X / by p ^ {y'^{p),^{p))i which satisfies y^{y,0) = y^. More- 
over, dy^{y,0) = dy'^{y,0), so g^^{y,0) = 0. The condition g^^{y,0) = 1 is not 
affected by changes in the y-coordinate, so defines the required tubular 

neighbourhood. ■ 
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